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ABSTRACT 


This  thesis  investigates  how  and  to  what  effectiveness  the  03b  commercial 
satellite  constellation  could  be  used  for  geolocation  of  Ka-band  EMI  sources  in  support 
of  the  DOD.  Review  of  commonly  used  geolocation  techniques  for  suitability  and 
comparison  of  those  with  the  03b  constellation  characteristics  shows  that  a  new  method 
of  geolocation  is  necessary  and  possible.  A  method  using  the  Doppler  effect  with 
frequency  data  from  a  single  antenna  was  then  created  that  is  compatible  with  03b.  This 
method  uses  the  received  frequency  of  the  jammer  over  time  to  detect  the  base  frequency, 
and  then  compares  the  received  frequency  to  that  of  simulated  emitters  at  known 
locations  in  order  to  provide  a  geolocation  for  EMI  emitters.  This  was  modeled  to  be 
accurate  within  16  km  throughout  03b ’s  service  area.  This  level  of  accuracy  would 
provide  the  ability  to  mitigate  the  interference  or  decrease  a  search  area  for  assets  with 
higher  capabilities  in  order  to  increase  their  efficiency  of  tasking/use.  This  method  is  of 
further  benefit  to  the  DOD  due  to  its  potential  to  be  low  cost,  be  maintained  as  an  organic 
capability  by  the  units,  and  decrease  the  time  necessary  to  reach  a  conclusion  when 
working  through  the  Joint  Spectrum  Interference  Resolution  process. 
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I. 


INTRODUCTION 


The  primary  research  question  for  this  thesis  was  how,  and  to  what  effectiveness, 
could  the  03b  commercial  satellite  constellation  be  used  for  geolocation  of  Ka-band  EMI 
sources  in  support  of  the  DOD?  There  has  been  an  increasing  move  of  SATCOM  users, 
including  military,  to  the  Ka-band  in  order  to  avoid  congestion  and  support  high  data 
rates.  The  military  also  already  relies  heavily  on  commercial  SATCOM  to  meet 
bandwidth  requirements.  Those  tendencies,  combined  with  the  need  for  geolocation 
services  in  order  to  preserve  the  SATCOM  capability,  and  the  possibility  that  such 
services  could  be  provided  cheaply  and  directly  to  the  user,  made  this  an  interesting 
research  focus  with  great  potential  benefit  to  the  DOD. 

In  breaking  the  primary  research  question  into  pieces  and  developing  the  subject, 
it  is  first  important  to  develop  what  SATCOM  EMI  is  and  why  it  is  an  issue.  These  will 
be  addressed  in  the  remainder  of  Chapter  I,  following  an  overview  of  satellite 
communications. 

Next,  the  case  for  why  geolocation  of  the  source  of  EMI  is  important  must  be 
made.  Geolocation  operations  in  the  Ka-band  present  particular  challenges  that  must  be 
addressed  as  well,  hollowing  those  issues,  an  overview  of  existing  geolocation  techniques 
must  be  made.  The  techniques’  suitability  for  use  will  then  be  determined.  These 
subsidiary  topics  will  be  covered  in  Chapter  II. 

A  technique  must  then  be  chosen  or  developed  and  applied  to  the  03b 
constellation  with  its  specific  traits  considered.  As  the  03b  constellation  is  unique  among 
communications  constellations  currently  in  use,  its  characteristics  will  be  discussed  in 
Chapter  III. 

The  development  or  selection  of  a  geolocation  technique  and  characterization  of 
the  accuracy  with  which  that  method  is  able  to  perform  geolocation  of  an  interference 
source  must  be  made.  Chapter  IV  of  this  thesis  will  do  that,  with  a  different  perspective 
on  how  to  use  Doppler  shift  as  determined  in  the  received  frequency  of  the  satellite.  This 
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is  evaluated  for  proof  of  eoneept  as  an  alternative  to  existing  methods  for  geoloeation  of 
EMI  emitters. 

Finally,  in  Chapter  V,  the  usefulness  of  the  suitable  teehnique(s)  for  geoloeation 
of  EMI  emitters  is  diseussed  with  the  emphasis  on  eapability  to  mitigate  the 
eommunieations’  interferenee  and  potential  benefits  to  the  DOD  if  03b  were  to 
implement  a  geoloeation  of  EMI  emitter  eapability. 

Chapter  VI  will  address  future  work  that  eould  be  done  to  improve  the  thesis  by 
inereasing  the  seope  and  utility,  improving  geoloeation  aeouraey,  or  investigating 
alternative  methods  to  address  the  researeh  questions. 

In  limiting  the  seope  of  the  thesis,  emphasis  was  plaeed  on  finding  a  teohnique(s) 
for  geoloeation  that  preferably  required  a  minimal  number  of  satellites  (ideally  one),  a 
minimum  number  of  antennas  (ideally  one),  and  is  a  “passive”  method — ^whieh  is  used  in 
this  ease  to  mean  it  is  only  reeeiving  signals.  This  thesis  will  be  best  understood  by 
readers  having  a  basie  understanding  of  satellite  eommunieations  (SATCOM)  and  orbital 
meehanies,  as  the  underlying  eoneepts  used  in  teehniques  for  geoloeation  of  EMI  emitters 
are  dependent  on  the  geometry  and/or  relative  motion  between  an  emitter  (the  EMI 
souree)  and  the  satellite  (EMI  reeeiver),  and  the  motion/eharaeteristies  of  radio  waves. 

A.  SATELLITE  COMMUNICATIONS  OVERVIEW 

The  idea  of  satellite  eommunieations  has  existed  sinee  at  least  1945,  when  Arthur 
C.  Clarke  was  published  diseussing  the  use  of  satellites  for  transmitting  television 
signals.!  The  idea  of  satellite  eommunieations  provided  a  solution  to  the  issue  of  our  non¬ 
flat  Earth  preventing  radio  waves  from  reaehing  distant  reeeivers.  Arthur  Clarke’s  idea 
was  greatly  expanded  upon  in  the  1950s  as  AT&T  took  an  interest  and  began  work  in 
earnest  at  making  the  idea  a  reality. 2  In  1957,  Sputnik  I  was  launehed,  not  speoifieally 
with  a  eommunieation  orientation,  but  as  the  first  man-made  satellite  an  important  leap 


!  David  G.  Whalen,  “Communications  Satellites  Short  History,”  NASA  History  Division,  November 
30,  2010,  http://history.nasa.gov/satcomhistory.html. 

2  Ibid. 
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forward  in  the  quest  for  satellite  eommunications.^  By  the  end  of  the  1950s,  NASA  had 
launched  communication-oriented  high-altitude  balloons,  which  served  as 
communications  reflectors,  and  the  U.S.  Department  of  Defense  (DOD)  was  researching 
developing  active  communications  satellites.^  In  the  early  1960s,  the  first  devoted 
communication  satellite,  TELSTAR,  was  launched.^  This  was  soon  followed  by  other 
programs  with  RELAY,  SYNCOM  and  EARLY  BIRD  all  by  the  middle  of  the  1960s.^ 
Since  then,  satellite  communications  have  grown  tremendously  and  now  exist  in  all  major 
orbital  regimes:  low  earth  orbit  (LEO),  medium  earth  orbit  (MEO),  geosynchronous  orbit 
(GEO),  and  highly  elliptical  orbit  (HEO).  Satellites  are  capable  of  serving  customers 
globally  and  play  a  huge  role  in  meeting  communication  needs  for  all,  from  individuals  to 
companies  to  militaries. 

Satellite  communications  are  conducted  by  transmission  and  reception  of  radio 
waves.  The  radio  waves  contain  information,  in  the  form  of  a  modulation  on  a  carrier 
wave.  This  is  transmitted  from  a  source  to  the  satellite,  which  receives  that  signal  and 
then  retransmits  it  to  the  intended  recipient.  Eor  our  look  at  satellite  communication,  we 
will  start  with  basic  satellite  communications  architecture.  The  architecture  can  be 
thought  of  in  three  segments.  These  are  ground,  space,  and  the  communications  link 
between  them.  The  “ground”  segment  more  accurately  refers  to  Earth-based  terminals, 
which  could  include  stationary  or  mobile  terminals  traversing  land,  sea,  or  air. 
“Terminal”  encompasses  the  necessary  equipment/hardware  and  software  required  to 
transmit  to  or  receive  communications.^  The  communication  link  portion  of  the 
architecture  joins  the  ground  and  space  piece  with  information  transmitted  via  radio 
waves,  and  includes  the  characteristics  of  the  environment  through  which  it  transits 
(distance,  rain,  atmosphere,  etc.).  This  link  includes  communication  of  both  “mission 
data”  (data  that  serves  the  purpose  the  satellite  was  designed  for  such  as  communications 

3  Ibid. 

Whalen,  “Communications  Satellites  Short  History.” 

5  Ibid. 

6  Ibid. 

^  Wiley  J.  Larson  and  James  R.  Wertz,  Space  Mission  Analysis  and  Design,  3rd.  edition  (Hawthorne, 
CA:  Microcosm  Press,  1999),  534. 
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or  imagery)  as  well  as  traeking,  telemetry,  and  eommand  (TT&C)  data  (for  the  purpose  of 
eontrolling  and  operating  the  satellite).^  The  term  for  information  traveling  from  ground 
to  spaee  is  “uplink,”  whereas  if  the  direetion  of  travel  is  from  spaee  to  ground,  the  term 
used  is  “downlink.”  Finally,  the  spaee  segment  eonsists  of  the  satellite  or  satellites, 
designed  to  reeeive  and  transmit  information  to/from  the  ground  segment  via  the 
eommunieation  link  segment.  Depending  on  the  satellite,  the  transmission  to  the  ground 
may  oecur  after  steps  that  eould  inelude  simple  amplifieation,  filtering,  and  translation  to 
a  downlink  frequency,  or  more  advanced  processing  to  include  error  correction  and  signal 
regeneration. 

To  understand  how  these  pieces  impact  communications,  the  link  Equation  (1.1) 
(commonly  referred  to  as  link  budget)  characterizes  the  strength  of  the  link,  typically  as 
the  signal-to-noise  ratio,  based  on  the  factors  that  affect  it.  Eh/No,  the  signal-to-noise  ratio 
normalized  per  bit,  directly  relates  to  the  probability  of  bit  error  through  modulation  and 
coding.9  On  a  simpler  level,  a  higher  Eh/No  means  fewer  errors  and  better  quality.  While 
the  link  equation  can  be  found  in  varying  forms,  one  that  is  useful  for  a  basic 
understanding  is  taken  from  Space  Mission  Analysis  and  Design  {SMAD)}^ 

EJK  =  tseefootnote]ll  (1.1) 

kT^R 

What  helps  increase  the  signal-to-noise  ratio  is  increased  transmitter  power  {P), 
transmit  antenna  gain  {Gt),  receiver  antenna  gain  (Gr),  or  reduced  free  space  loss  (T^), 
transmission  path  loss  {La),  line  loss  (T/),  data  rate  {R),  and  system  noise  temperature 
(rs)!^.  The  {k)  in  the  equation  is  Boltzmann’s  constant.  Not  all  of  these  items  are 
adjustable.  Free  space  loss  is  a  function  of  the  distance  from  the  transmitter  to  the 
receiver,  and  therefore  is  not  easily  changed  when  dealing  with  satellite  communications. 
The  Earth’s  atmosphere,  rain,  foliage,  and  similar  obstructions  result  in  transmission  path 

^  Larson  and  Wertz,  Space  Mission  Analysis  and  Design,  534. 

9  Ibid.,  561. 

10  Ibid.,  551. 

11  Ibid.,  551. 

12  Ibid.,  551. 
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loss,  particularly  at  higher  frequencies  of  communication  such  as  the  Ka-band,  which  is 
where  this  thesis  will  be  concerned.  Line  loss  is  a  function  of  the  hardware  used  in  the 
ground  and  space  segments  of  the  communication  system,  and  is  typically  extremely 
small  in  comparison  to  free  space  loss.  The  system  noise  temperature  is  an  equivalent 
temperature  representing  the  summation  of  the  degrading  effect  resulting  from  various 
sources  of  environmental  noise  both  natural  (e.g.,  sun  or  Earth’s  presence  in  antenna 
beam)  and  man-made  (buildings).  1 3  Gain  for  any  antenna  (transmitter  or  receiver)  is  a 
function  of  antenna  size,  wavelength,  and  efficiency.  While  the  former  factors  are  often 
difficult  to  change  as  characteristics  of  the  environment,  transmitter,  receiver;  transmitter 
power  sometimes  can  be  changed  or  increased  within  system  constraints  to  provide  an 
increased  Eb/No.  Data  rate  may  also  be  decreased  to  provide  increased  Eb/No  and 
improved  communications  quality. 

B.  EMI  AND  JAMMING  OVERVIEW 

In  satellite  communications  today  a  common  problem  for  both  civilian  and 
military  users  is  electromagnetic  interference  (EMI).  The  Chairman  of  the  Joint  Chiefs  of 
Staff  Manual  (CJCSM)  3320.02D  -  Joint  Spectrum  Interference  Resolution  (JSIR) 
Procedures  defines  EMI  as  “any  electromagnetic  (EM)  disturbance  that  interrupts, 
obstructs,  or  otherwise  degrades  or  limits  the  effective  performance  of  electronics  or 
electrical  equipment.”i4  The  manual  goes  on  to  state  that  the  causes  of  EMI  may  be  either 
unintentional  or  intentionally  This  should  be  added  to  the  definition  of  interference  by 
Carlson  and  Crilly,  which  clarifies  interference  is  “contamination  by  extraneous  signals 
from  human  sources,”!^  regardless  of  whether  it  is  intentional  or  not.  Other  causes  of 
degradation  include  noise,  which  has  been  defined  as  “random  and  unpredictable 
electrical  signals  produced  by  natural  processes  both  internal  and  external  to  the  system” 
and  distortion,  defined  as  “waveform  perturbation  caused  by  imperfect  response  of  the 

13  Larson  and  Wertz,  Space  Mission  Analysis  and  Design,  556. 

14  “Joint  Spectrum  Interference  Resolution  (JSIR)  Procedures,”  Chairman  of  the  Joint  Chiefs  of  Staff 
Manual  (CJCSM  3320.02D),  June  3,  2013,  A-1. 

ly  Ibid.,  A-1. 

1^  A.  Bruce  Carlson  and  Paul  B.  Crilly,  Communication  Systems:  An  Introduction  to  Signals  and  Noise 
in  Electrical  Communications,  5*  edition  (New  York:  McGraw-Hill,  2011),  4. 


5 


system  to  the  desired  signal  itself.’’^^  When  undesired  radio  waves/energy  are  reeeived  by 
the  satellite  the  desired  signal  and  information  from  the  undesired  signal  and  elutter 
beeome  harder  to  distinguish.  If  the  undesired  and  desired  signals  are  unable  to  be 
eompletely  and  aeeurately  separated  from  eaeh  other,  the  quality  (mateh  of  what  is 
reeeived  to  what  was  sent)  of  the  transmission  is  negatively  impaeted.  The  seale  of  the 
impaet  depends  largely  on  the  relative  strength  of  the  various  signals  involved  and  the 
ability  of  processing,  modulations,  and  coding/error  correction  methods  available  to  help 
distinguish  desired  from  undesired  signals. 

Since  out  of  the  three  degradation  sources  (distortion,  noise,  and  interference) 
interference  is  the  only  one  caused  by  humans,  that  will  be  our  focus — the  others  can  be 
thought  of  as  the  cost  of  doing  business  in  satellite  communications  and  can  also  be 
corrected  to  some  extent.  Causes  of  interference  are  numerous.  Incorrect  frequency  use 
by  communications  users,  other  transmitters  in  the  antenna  beam  (ranging  from  friendly 
to  adversarial),  equipment  improperly  set  up  or  in  proximity  to  other  electrical 
equipment/electronics,  or  improperly  pointed  antennas,  along  with  anything  else  in  the 
form  of  man-made  items  or  use  of  them  either  properly  or  not.  The  term  jamming  in  this 
thesis  will  be  defined  as  EMI  caused  intentionally  by  a  hostile  emitter. 

1.  Why  EMI/Jamming  Is  an  Issue  for  the  DOD 

Satellite  communications  are  extremely  important  to  the  Department  of  Defense. 
A  report  by  the  Government  Accountability  Office  has  stated  that  “DOD’s  use  of 
SATCOM  is  critical  to  military  operations  worldwide  and  its  dependency  is  expected  to 
increase  over  the  next  decade.”!^  The  same  report  stated  that  between  2000  and  2011  the 
DOD  increased  commercial  SATCOM  use  by  800%;  presumably  reflecting  close  to  the 
total  increase  in  SATCOM  use  by  the  DOD  due  to  the  small  percentage  of  SATCOM  met 
by  DOD-owned  assets.  The  inherent  nature  of  the  United  States  military  as  a 

1  ^  Carlson  and  Crilly,  Communication  Systems:  An  Introduction,  4. 

Government  Accountability  Office.  Defense  Satellite  Communications:  DOD  Needs  Additional 
Information  to  Improve  Procurements  (GAO-15-459)  (Washington,  DC:  Government  Accountability 
Office,  2015),  19. 

19  Ibid.,  9. 
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deployable,  expeditionary  foree — expected  to  execute  missions  at  locations  around  the 
globe,  often  in  remote  or  undeveloped  locations — results  in  a  need  for  a  communications 
structure  that  can  accommodate  those  missions  on  minimal  notice.  Satellite 
communications  contribute  to  the  capability  to  meet  that  need.  With  no  permanent 
ground  infrastructure  required,  small  terminals  can  be  deployed  and  set  up  or  even  carried 
by  ship,  vehicle,  or  personnel.  These  allow  command  and  control,  communications,  and 
intelligence  to  occur  unimpeded  and  be  distributed  in  an  expedient  manner  anywhere  in 
the  world.  The  need  for  satellite  communications  capacity  has  drastically  risen  with  the 
increased  demand  for  intelligence,  surveillance  and  reconnaissance  (ISR)  products  such 
as  high  resolution  imagery  and  video  used  for  battle  space  awareness  and  targeting.  Any 
reduction  in  availability  of  satellite  communications  due  to  jamming  could  slow  or  even 
severely  degrade  military  operations. 

2.  How  Jamming  Occurs 

In  order  to  understand  how  jamming  typically  occurs  some  geometry  must  be 
considered.  Before  that,  antenna  gain  as  it  relates  to  beamwidth  must  be  addressed. 
Antenna  gain  at  its  most  basic  is  a  numerical  description  of  the  ability  to  concentrate 
signal  energy  that  would  otherwise  be  distributed  in  all  directions.  For  a  common 
parabolic  antenna  with  a  circular  aperture,  an  antenna  with  a  larger  diameter  or  operating 
at  high  frequencies  (shorter  wavelengths)  is  able  to  concentrate  energy  better  than  a 
smaller  one  or  one  operating  at  lower  frequencies  (longer  wavelengths).  The  projection  of 
this  area  of  concentrated  signal  strength  onto  the  Earth  (for  a  satellite-based  antenna)  can 
be  called  the  beam  coverage  area  (BCA)  or  “footprint.”  If  the  antenna  in  question  is 
located  on  the  Earth  and  aimed  at  a  satellite,  then  the  BCA  would  fall  in  space.  The  gain 
is  highest  in  the  center  of  this  area  (called  the  boresight)  and  then  falls  off  in  a  lobe 
pattern  at  increasing  distances  from  center.  Due  to  this  pattern,  the  width  of  the  beam 
must  also  specify  the  level  of  gain  used  to  characterize  that  width.  The  -3  dB  point  from 
peak  gain  is  typically  used,  and  the  diameter  (measured  in  degrees  to  accommodate 
satellite  altitudes)  is  aptly  named  the  “half-power  beamwidth.”  The  beam  pattern  for  an 
03b  payload  antenna  demonstrates  this  concept  in  Eigure  1 . 
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Figure  1 .  03b  Satellite  Payload  Antenna  Pattern  with  Half-Power  Beamwidth 

and  Side  Lobes.^o 


— ^RxtheU  cut 
— ^Rxphi  cot 


Moving  to  the  geometry  portion  of  jamming,  in  order  for  effective  reception  of  a 
signal  to  take  place,  a  transmitter  should  be  located  within  the  BCA  at  a  location  with 
suitable  gain  to  meet  the  objective.  This  is  true  regardless  of  whether  the  transmitter  is 
located  at  the  satellite  or  ground  terminal,  and  regardless  of  whether  the  transmitter  is  the 
intended  one  or  a  hostile  emitter — the  only  difference  is  that  the  objective  shifts  from 
closing  a  link  to  preventing  closure  of  a  link  if  the  transmitter  is  hostile.  With  SATCOM, 
the  reason  the  uplink  is  typically  jammed  instead  of  the  downlink  (a  notable  exception 
being  GPS  jamming,  which  is  not  covered  in  this  thesis)  is  that  to  jam  the  uplink  the 
transmitter  can  be  on  the  earth,  and  to  jam  the  downlink  it  would  have  to  be  somewhere 
above  the  receiver  in  the  beam — such  as  space. 

While  satellite  communications  offer  significant  advantages  over  some  other 
forms  of  communication,  they  still  have  vulnerabilities,  many  resulting  from  the  same 
things  that  make  them  so  useful.  A  fundamental  disadvantage  of  satellite  communications 
is  that  a  critical  part,  the  satellite,  has  to  be  in  space.  This  means  that  it  had  to  be  placed 
in  orbit,  and  that  causes  constraints  on  weight,  size,  and  power  available.  Antenna 


Adapted  from  J.  J.  Shaw,  “03b  Networks  and  Why  Not?”  (presentation,  Naval  Postgraduate  School, 
Monterey,  May  7,  2015). 
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diameter  (size)  is  an  important  factor  in  the  link  budget  (1.1)  that  significantly  impacts 
the  signal-to-noise  ratio.  Satellites  are  also  expensive.  Furthermore,  the  purpose  of 
SATCOM  as  the  DOD  uses  it  is  often  to  provide  communications  to  a  user  who  is  using  a 
small  terminal  in  order  to  maintain  deployable,  expeditionary,  and  mobile 
communications.  The  issue  becomes  that  a  jammer  on  the  ground  has  a  dramatic 
advantage  over  the  satellite  when  considering  cost  since  there  is  no  cost  per  pound  to  lift 
it  to  space,  and  can  be  comparable  to  or  surpass  the  friendly  force  ground  transmitter  in 
important  aspects  affecting  gain.  The  jammer  can  easily  have  a  large  antenna  and  power 
supply  as  even  a  mobile  jammer  is  only  limited  by  the  size  and/or  number  of  vehicles 
required  to  carry  it.  A  stationary  jammer  faces  even  less  restriction.  While  Ka-band 
technology  is  more  complex  than  lower  frequencies,  such  as  UHF,  the  technology  is 
widespread  and  available  for  purchase  in  C  and  Ku  bands  (an  internet  search  for  “Ku- 
band  SATCOM  jammer”  can  demonstrate  this)  and  could  reasonably  be  expected  to 
progress  to  Ka-band  as  the  band  becomes  more  commonly  used,  and  therefore  a  larger 
target. 21  In  countries  where  Ka-band  coverage  exists  a  television  news  van  or  other 
terminal  (whether  purchased  or  legally  or  illegally  acquired)  set  up  to  transmit  by  satellite 
could  easily  be  repurposed  as  a  jammer  by  the  adversary,  or  inadvertently  become  one 
through  user  error  by  news  crews  covering  events  in  the  area  or  due  to  equipment 
malfunction.22  Numerous  websites  show  entire  fleets  of  Ka-band  operating  news  vans, 
many  of  which  are  advertised  as  low-cost.23 

While  jamming  technology  is  available  and  a  relatively  low  cost  method  of 
interfering  with  satellite  communications,  to  be  most  effective  it  still  needs  to  be  located 
within  the  BCA  of  the  satellite’s  antenna,  and  have  the  transmit  antenna  oriented  at  the 
receiver  satellite.  Even  for  Ka-band  frequency  satellites  with  relatively  small  spot  beams 


21  Ronald  C.  Wilgenbusch  and  Alan  Heisig,  “Command  and  Control  Vulnerabilities  to 
Communications  Jamming,”  Joint  Forces  Quarterly,  no.  69  (April  2013):  58. 

22  Ibid.,  61. 

23  Skylogic.  NEWSSPOTTER.  (accessed  May  17,  2016). 
http://www.skylogic.it/?page_id=2309&lang=en;  C-Com  Satellite  Systems  Inc.  Ka-Band,  (accessed  May 
17,  2016).  http://www.c-comsat.com/solutions/ka-band/;  Frontline  Communications.  Mobile  Ka-Band 
Solution  Debuts  for  Satellite  Newsgathering.  August  21,  2013. 
http://www.fronthnecomm.com/news/KaWriteUp.pdf 
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there  ean  be  signifieant  antenna  gain  over  hundreds  of  km,  with  an  example  being  03b, 
whieh  provides  greater  than  25  dB  gain  for  an  area  more  than  700  km  wide.^^  Without 
publishing  where  these  beams  are  pointed  or  unit  loeations,  if  hostile  forees  happen  to 
know  where  DOD  forces  are  located/operating  within  hundreds  of  km,  which  they  likely 
will  if  the  operations  concern  them,  then  they  can  easily  fall  within  the  BCA.  Another 
scenario  would  be  that  the  actor  is  detecting  down-link  signals  and  could  then  assume 
they  are  within  the  BCA.  Neither  scenario  is  unreasonable.  For  example,  at  a  500  km 
diameter,  a  beam  would  cover  the  majority  of  Syria.  A  1000  km  wide  beam,  on  the  other 
hand,  would  cover  nearly  all  of  Iraq.  Finding  the  satellites  that  are  likely  used  for 
communications  are  also  relatively  simple,  as  DOD  communications  satellites  such  as 
Wideband  Global  SATCOM  (WGS)  are  not  classified  and  there  is  a  wide  number  of 
amateur  astronomy  websites  available  which  post  ephemeris  data  for  both  commercial 
and  military  satellites.  Furthermore,  due  to  satellite  communications  being  restricted  to 
line  of  sight  and  the  limited  number  of  Ka-band  communications  satellites,  there  would 
be  a  small  number  of  options  for  which  satellites  could  be  in  use  for  a  given  area.  This 
would  then  allow  the  hostile  force  to  orient  the  jammer  antenna  on  the  communications 
satellite  to  optimize  their  chances  of  successful  interference. 

C.  CONCLUSION 

With  their  use  on  a  global  scale  by  the  military  for  important  functions,  satellite 
communications  are  a  likely  target  for  adversaries.^^  In  the  past  years,  intentional 
jamming  incidents  have  risen  rapidly,  with  numerous  incidents  of  government  or  political 
bodies  in  the  Middle  East  intentionally  jamming  various  satellites  having  been  widely 
reported.26  One  article  citing  a  high-level  Eutelsat  employee  placed  the  number  of 
intentional  interference  incidents  Eutelsat  experienced  at  54  in  2010,  twice  that  in  2011, 


24  J.  J.  Shaw,  “03b  Networks  and  Why  Not?” 

25  Wilgenbusch  and  Heisig,  “Command  and  Control  Vulnerabilities  to  Communieations  Jamming,” 
56-63. 

26  Ibid.,  58. 
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and  three  times  the  2011  figure  in  11  months  of  2012.27  This  growing  threat  of 
eommunieations  jamming  is  one  that  will  need  to  be  addressed;  while  geoloeation  of  the 
jammer  does  not  prevent  jamming  from  oeeurring,  it  ean  help  resolve  the  threat  onee  it 
oeeurs. 


27  Anne-Wainscott  Sargent,  “Fighting  Satellite  Interference  on  All  Fronts,”  Via  Satellite  - 
SatelliteToday.com,  March  1,  2013,  http://www.satellitetoday.com/pubhcations/via-satellite- 
magazine/features/20 1 3/03/0  l/fighting-satellite-interference-on-all-fronts/. 
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II.  GEOLOCATION  OF  EMI  SOURCES 


Chapter  I  gave  a  quick  overview  of  satellite  communications,  EMI,  and  why  this 
is  a  problem  for  the  DOD.  Chapter  II  will  first  address  the  importance  of  the  ability  to 
locate  the  source  of  EMI.  It  will  next  cover  some  specific  characteristics  of  Ka-band, 
which  are  important  to  the  focus  of  this  particular  thesis.  Einally,  it  will  discuss  some 
existing  methods  of  geolocation  and  why  they  are  not  ideal  for  use  with  Ka-band 
communications. 

A.  WHY  FINDING  EMI  SOURCES  IS  IMPORTANT 

In  CJCSM  3320.02D,  geolocation/direction  finding  is  listed  as  a  part  of  the  JSIR 
process  following  initial  detection,  verification,  characterization,  and  reporting. 
Geolocation  is  important  because  the  resolution  process  depends  on  the  location, 
ownership,  and  users  of  the  problem  emitter  and  the  impacted  receiver.28  Ownership  and 
users  would  be  difficult  to  identify  without  knowing  the  location  of  the  emitter.  Also, 
although  geolocation  occurs  after  characterization  in  the  process,  knowing  a  location 
could  confirm  an  interference  incident  was  hostile,  neutral,  or  friendly.  Einally,  even  if 
geolocation  fails  to  provide  insight  as  to  the  ownership  or  users  of  the  problem  emitter,  a 
location  of  the  emitter  can  enable  mitigation  of  the  interference  and  restoration  of  the 
degraded  resource — topics  which  will  be  addressed  in  greater  depth  later.  When  dealing 
with  interference  incidents  outside  of  simulations  such  as  in  this  thesis  it  is  important  to 
ensure  all  information  pertaining  to  the  incident  and  any  associated  geolocation  is 
classified  to  the  appropriate  level — guidance  may  be  found  in  CJCSM  3320.02D 
Enclosure  D.29 

B.  CURRENT  GEOLOCATION  TECHNIQUES  OVERVIEW 

In  the  search  for  techniques  that  could  be  utilized  by  03b  communications 
satellites,  the  first  step  was  to  learn  about  existing  geolocation  techniques  and  evaluate  if 

“Joint  Spectrum  Interference  Resolution  (JSIR)  Procedures,”  A-7-A-13. 

29  Ibid.,  A-6,  A-7. 


13 


they  were  suitable  for  use.  Given  the  variety  of  methods,  these  ean  be  organized  by 
similarity  in  requirements,  sueh  as  number  of  satellites/reeeivers. 

1.  Multiple  Satellite  Methods 

Two  of  the  most  common  methods  of  geolocation  are  time  difference  of  arrival 
(TDOA)  and  frequency  difference  of  arrival  (FDOA),  which  both  require  multiple 
satellites.  The  required  geometry  for  both  methods  is  the  same,  and  consists  of  two 
satellites  receiving  the  interference  signal  simultaneously.  With  TDOA  a  specific  pulse  or 
portion  of  the  signal  must  be  able  to  be  uniquely  identified  in  order  that  a  time  difference 
in  when  that  pulse  is  received  by  the  two  different  satellites  may  be  measured.  The  time 
difference  then  correlates  to  a  difference  in  distances  between  the  respective  satellites  and 
transmitter.  A  key  aspect  is  that  while  the  difference  in  distances  is  known,  the 
magnitudes  of  each  distance  is  not — as  the  receiver  does  not  know  when  the  signal  was 
originally  transmitted.  This  difference  in  time/distances  is  constant  over  a  number  of 
possible  locations  which  form  a  line  of  position  (LOP),  known  as  an  isochrone,  for 
“same-time. ”30  A  single  line  of  position  is  inadequate  for  a  geolocation,  but  when  various 
TDOAs  are  collected  or  another  method  is  used  to  obtain  an  LOP  then  positional 
ambiguity  can  be  resolved  and  accuracy  is  increased.  This  method  is  well  discussed  in 
Dr.  Herschel  Loomis’  paper  Geolocation  of  Electromagnetic  Emitters?^  The  paper  also 
discusses  a  similar  TDOA  approach,  the  “Interferometric  Approach,”  in  which  if  the 
distance  between  satellites  is  small  compared  to  the  distance  to  the  transmitter,  the  math 
may  be  simplified  and  results  in  an  angle  giving  the  isochrone/LOP.32  The  resultant  line 
is  noted  to  be  the  “same  line  which  defined  the  asymptotes  of  the  hyperbolic 
isochrones. ”33  An  example  of  TDOA  isochrones  from  Geolocation  of  Electromagnetic 


30  Herschel  H.  Loomis,  Geolocation  of  Electromagnetic  Emitters  (NPS-EC-00-003),  (Monterey: 
Naval  Postgraduate  School,  Revised  October  2007),  1 . 

31  Ibid. 

32  Ibid.,  6. 

33  Ibid.,  5-6. 
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Emitters  is  shown  in  Figure  2,  where  the  satellite  would  be  moving  left  to  right  along  the 

x-axis.34 


Figure  2.  Family  of  Hyperbolic  TDOA  Isochrons. ^5 


FDOA  is  often  used  as  the  complement  to  the  TDOA  LOP.  The  reason  that  these 
methods  work  well  together  is  that  while  TDOA  provides  an  LOP  that  is  generally  an 
angle  along  which  the  emitter  lies,  FDOA  provides  a  range.36  The  FDOA/range  LOP 
intersects  the  TDOA  LOP  to  provide  a  geolocation.  FDOA  once  again  requires  two 
satellites  that  receive  the  signal  from  the  transmitter  simultaneously.  FDOA  also  requires 
that  the  satellites  are  moving  at  sufficient  speed  to  generate  a  Doppler  shift  significant  in 
comparison  to  the  ability  to  measure  frequency  (so  it  is  not  lost  in  measurement  error). 
Doppler  shift  {Ef,  Hz)  as  shown  in  Equation  (1.2)  is  a  result  of  the  velocity  (v,  in  m/s) 

Loomis,  Geolocation  of  Electromagnetic  Emitters,  5. 

Source:  Loomis,  Geolocation  of  Electromagnetic  Emitters,  5. 

36  Darko  Musicki  and  Wolfgang  Koch.  “Geolocation  using  TDOA  and  FDOA  Measurements,”  in 
Proceedings  of  the  11th  International  Conference  on  Information  Fusion,  (Cologne:  IEEE,  2008),  1987- 
1994. 
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component  of  the  reeeiver  along  veetor  r  eonneeting  the  transmitter  and  reeeiver.  This  is 
then  multiplied  by  the  frequeney  (Hz)  and  divided  by  the  speed  of  light  (c,  m/s). 
Convention  is  that  if  the  range  from  reeeiver  to  transmitter  is  deereasing  v  and  A/^  are 
positive. 

(1.2) 

c  r  c  c 


This  is  beeause  if  the  transmitter  and  reeeiver  are  elosing  in  distanee  then  signal  waves 
are  eompressed  whieh  results  in  a  frequeney  upshift.  If  the  transmitter  and  reeeiver  are 
inereasing  in  respeetive  range  then  frequeney  downshift  oeeurs  as  the  signal  waves  are 
spread  more.  If  there  was  no  relative  motion  the  Doppler  shift  would  be  0.  The  FDOA 
Equation  (1.3),  as  taken  from  Geolocation  of  Electromagnetic  Emitters,  for  frequeneies  // 
and  f2  reeeived  by  reeeivers  1  and  2,  respeetively,  demonstrates  in  two  dimensions  how 
the  known  reeeiver  veloeities  and  measured  FDOA  would  relate  to  a  LOP  eharaeterized 
by  X  and  y  eoordinates.^^  This  equation  represents  a  situation  where  the  reeeivers  are 
loeated  on  the  x-axis  at  {s,0)  and  {-s,0),  eaeh  moving  at  veloeity  {vx,Vy),  and  the  emitter  is 
at  {xe,  ye).38 


¥E>OA=f-f 


/ 

c 


^x{^e+s)  +  ^yye 
+  ye' 


A{^e-^)  +  fye 

^l{x^-sf+yf  , 


[see  footnote]^^ 


(1.3) 


This  equation  is  the  result  of  taking  the  differenee  of  the  Doppler  shift  reeeived  by 
eaeh  reeeiver  (as  seen  in  Equation  (1.2)  for  a  single  reeeiver),  and  then  expanding  the 
veloeity  of  the  reeeivers  and  r  veetors  between  eaeh  reeeiver  and  emitter  into 
eomponents.  The  solution  to  Equation  (1.3)  results  in  an  LOP  giving  range.  When  used  in 
eonjunetion  with  TDOA  giving  bearing,  or  other  LOPs  the  ambiguity  is  resolved  and 
results  in  a  geoloeation.  This  is  demonstrated  in  Figure  3  taken  from  Geolocation  of 


Loomis,  Geolocation  of  Electromagnetic  Emitter,  6-7. 

38  Ibid.,  7. 

39  Ibid.,  7. 
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Electromagnetic  Emitters,  which  shows  how  FDOA  contours  (blue)  complement  the 
TDOA  isoehrones  (red)  to  create  interseetions  that  would  identify  a  position.^o 


Figures.  TDOA  and  FDOA  Contours. 


2.  Multiple  Antenna  Methods 

For  methods  of  geoloeation  using  multiple  antennas  a  eommon  method  is  angle  of 
arrival  (AoA),  basically  using  the  interferometric  approach  to  TDOA  as  discussed 
previously,  but  with  the  reeeivers  as  multiple  antennas  on  a  single  satellite  as  opposed  to 
on  multiple  satellites.  An  alternative  method  whieh  has  been  modeled  and  may  eurrently 
be  undergoing  real  world  testing  is  using  steerable  beams  to  ereate  a  pattern  of  footprints 
from  which  signal  strength  differenees  ereated  by  the  emitter’s  varying  loeation  in  each 
beam  can  be  used  to  ealeulate  lines  of  position.42 


“^0  Loomis,  Geolocation  of  Electromagnetic  Emitters,  10. 

Source:  Loomis,  Geolocation  of  Electromagnetic  Emitters,  10. 

Brian  C.  Fredrick,  “Geoloeation  of  Source  Interference  From  A  Single  Satellite  With  Multiple 
Antennas”  (master’s  thesis.  Naval  Postgraduate  School,  2014). 
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3.  Single  Antenna  Methods 

The  primary  method  of  using  a  single  antenna  for  geoloeation  of  an  emitter  is 
interferometry.  This  method  is  analogous  to  TDOA/AoA,  as  phase  difference  is  directly 
related  to  time  difference.  The  difference  is  simply  in  what  property  is  able  to  be 
measured  based  on  the  situation.  For  single-antenna  interferometry  typically  a  phased 
array  antenna  is  used  which  has  elements  spaced  extremely  close  together,  at  less  than 
half  the  wavelength  of  the  transmitted  frequency.  The  spacing  is  important  as  there  would 
otherwise  be  ambiguity  as  to  whether  the  measured  phase  difference  was  correct  or  off  by 
a  multiple  of  2n.  By  measurement  of  the  phase  difference  in  the  signal  as  received  by  two 
elements  on  the  array,  and  using  the  known  wavelength,  an  angle  may  be  computed.  A 
phase  difference  that  correlates  to  the  spacing  between  elements  would  mean  that  the  two 
elements  formed  a  line  pointing  towards  the  transmitter,  whereas  no  phase  difference 
would  mean  that  the  elements  lie  on  a  line  perpendicular  to  a  line  to  the  transmitter  (a 
bearing  line).  Of  course,  in  between  these  boundary  scenarios  lie  phase  differences 
correlating  to  a  number  of  other  transmitter  bearings.  The  bearing  line  calculated,  as  with 
the  interferometric  approach  to  TDOA,  is  the  asymptote  of  a  hyperbola  with  the  foci  at 
one  of  the  receiver  elements.  Which  asymptote  the  transmitter  lies  on  is  initially 
ambiguous  as  the  sign  of  the  angle  calculated  is  unknown,  until  the  method  is  repeated 
with  other  elements.  This  method  and  the  associated  math  developing  the  asymptote  lines 
for  the  hyperbolas  are  detailed  in  Jiadio  Interferometric  Angle  of  Arrival  Estimations^ 

An  alternative  method  of  geoloeation  using  a  single  antenna  is  through  the  use  of 
the  Doppler  effect.  The  significant  advantage  of  the  use  of  the  Doppler  effect  is  that  it 
requires  only  a  single  antenna  on  a  single  satellite,  and  does  not  require  a  priori 
transmitter  frequency  knowledge.  The  Doppler  effect  varies  in  accordance  with  relative 
velocity  of  the  transmitter  to  receiver,  as  Equation  (1.2)  previously  showed.  The  relative 
velocity,  when  the  transmitter  is  stationary  and  the  orbit  is  known,  can  be  correlated  to 
physical  location. 


Isaac  Amundson  et  al.,  “Radio  Interferometric  Angle  of  Arrival  Estimation,”  in  Proceedings  of 
2010  European  Conference  on  Wireless  Sensor  Networks,  (Coimbra:  Springer,  2010). 
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C.  WHY  GEOLOCATING  EMI  SOURCES  IN  KA-BAND  IS  CHALLENGING 

The  military  today  requires  high  quantities  of  data  to  aecommodate  everything 
from  everyday  communieations  for  operations,  to  morale  boosting  applications,  to  large 
quantities  of  intelligence.  The  required  bandwidth  is  only  increasing,  as  discussed  in 
Chapter  I,  Section  B.l.  Other  frequency  bands,  such  as  C  and  Ku-band  generally  support 
lower  data  rates,  have  wider  beam  coverage  areas,  many  more  users,  all  resulting  in  them 
becoming  increasingly  congested.  This  is  going  to  result  in  the  military  increasingly 
shifting  to  Ka-band  to  facilitate  its  needs,  and  with  the  increased  shift  to  Ka-band  will 
come  the  need  for  geolocation  services  in  the  Ka-band. 

There  are  numerous  challenges  to  geolocating  Ka-band  EMI  sources.  Many  of  the 
same  characteristics  that  make  Ka-band  satellites’  use  beneficial,  such  as  their  difficulty 
to  jam  when  compared  with  lower  frequencies,  also  make  using  them  for  geolocation 
challenging.  Military  Ka-band  operates  with  typical  uplink  frequencies  of  about  30  GHz, 
as  seen  in  the  WGS  constellation.44  Commercial  Ka-band  is  close  to  that,  though 
typically  slightly  lower.  Use  of  higher  frequencies  provides  increased  bandwidth  which 
can  support  either  higher  levels  of  protection  through  frequency  spreading  or  hopping,  or 
support  much  higher  data  rates  than  lower  frequencies.  These  characteristics  and  others 
related  to  them  are  useful  in  making  communications  more  robust,  but  these  frequencies 
and  the  current  constellations  also  eliminate  several  existing  techniques  as  options  for 
geolocation  of  EMI  sources. 

I.  Number  of  Satellites 

Currently,  there  are  significantly  fewer  satellites  operating  in  Ka-band  than  in 
other  frequencies  such  as  C  or  Ku-band.  While  this  helps  with  congestion  and  reduces  the 
likelihood  of  unintentional  interference,  there  are  types  of  geolocation  methods  such  as 
EDO  A  and  TDOA  which  utilize  multiple  satellites.  Eor  these  methods,  each  satellite  must 
have  the  jammer  within  the  antenna’s  footprint  to  compare  information  about  the  signal 
as  it  is  received  by  other  satellites  at  the  same  time.  To  remedy  this  problem,  a  method 

Richard  A.  William,  and  Heywood  I.  Paul,  “Potential  Uses  of  the  Military  Ka-band  for  Wideband 
MILSATCOM  Systems,”  m  Military  Communications  Conference,  (Boston:  IEEE,  1998),  30-34. 
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able  to  provide  Ka-band  frequeney  geoloeation  using  a  single  satellite  and  single  antenna 
on  that  satellite  should  be  used  if  possible. 


2.  Antenna  Gain  and  Beamwidth 


The  higher  frequencies  that  fall  under  Ka-band  give  greater  antenna  gains  at 
smaller  sizes,  making  Ka-band  more  suitable  for  mobile/on-the-move  systems  requiring 
high  data  rates  than  lower  frequencies.  The  antenna  gain  Equation  (1.4)  taken  from 
(SMAD)  is  shown  below,  where  D  is  antenna  diameter  in  meters,  rj  is  an  efficiency  factor, 
and  A  is  wavelength  in  meters.  ^5 


G  = 


[see  footnote]^® 


(1.4) 


A  quick  comparison  of  the  30  GHz  Ka-band  uplink  to  a  12  GHz  Ku-band  signal 
with  a  1  meter  antenna  (for  ease  of  math)  using  the  above  equations  shows  the  Ka-band 
signal  has  a  gain  of  more  than  six  times  (about  8  dB)  greater  than  at  the  Ku-band 
frequency.  This  larger  gain  can  help  overcome  interference  or  otherwise  improve  the  link 
budget  to  better  meet  the  needs  of  the  user. 

The  higher  frequency  of  Ka-band  also  gives  a  smaller  beamwidth  for  a  given 
antenna  size.  The  half-power  (defined  as  a  3  dB  drop  in  signal  strength  from  beam  center) 
beamwidth  Equation  (1.5)  as  taken  from  SMAD  is  shown,  where  beamwidth  is  in  degrees, 
D  is  antenna  diameter  in  meters,  and /ghz  is  the  frequency  in  GHz. 


e  = 


[see  footnote] 


(1.5) 


Another  quick  comparison  of  a  30  GHz  signal  for  a  Ka-band  uplink  to  a  12  GHz 
Ku-band  signal  with  a  1  meter  antenna  shows  the  Ka-band  signal  has  a  half-power 


45  Larson  and  Wertz,  Space  Mission  Analysis  and  Design,  553. 

46  Ibid. 

47  Ibid.,  555. 

48  Ibid. 
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beamwidth  of  .7  degrees,  where  the  Ku-band  frequeney  has  a  half-power  beamwidth  of 
1.75  degrees.  The  differenee  may  seem  insignilieant,  but  at  GEO  one  degree  differenee  is 
over  600  km.  If  the  Ka-band  system  is  at  MEO  instead,  the  differenee  may  be  almost 
1000  km.  The  wider  beams  that  are  resultant  from  lower  frequeneies  (sueh  as  Ku-band) 
mean  that  a  ground-based  emitter’s  beam  is  more  likely  to  eneounter  multiple  GEO 
slots/satellites,  and  satellites  near  eaeh  other  in  GEO  are  more  likely  to  have  overlapping 
BCAs  on  the  ground.  Due  to  these  eonditions,  the  probability  of  multiple  satellites 
reeeiving  interferenee  from  a  eommon  transmitter  is  higher,  and  with  that,  the  ability  to 
eonduet  multiple  satellite  geoloeation  teehniques  is  also  more  likely.  The  same  eonditions 
and  opportunities  are  unlikely  for  a  Ka-band  system  due  to  the  mueh  smaller  beamwidths 
and  possibility  of  multiple  satellites  reeeiving  interferenee  from  a  eommon  transmitter. 

A  eommon  term  used  to  deseribe  the  small  and  typieally  steerable  beams  now 
frequently  used  by  Ka-band  and  high  throughput  satellites  (HTS)  is  “spot  beam,”  whieh 
is  in  eontrast  to  very  wide  beams  designed  to  eover  a  larger  area  of  the  Earth  (sometimes 
everything  in  the  line  of  sight)  without  steering.  This  smaller  beam  is  a  desired  and 
purposefully  engineered  aspeet  of  Ka-band  satellites,  as  it  provides  numerous  benefits.  A 
small,  steerable  beam  is  able  to  eoneentrate  gain  where  it  is  most  needed  (or  paid  for) 
while  wasting  less  on  areas  that  the  user  does  not  serve.  The  spot  beam  also  makes 
jamming  more  diffieult  due  to  the  fact  that  the  emitter  must  be  within  a  smaller 
geographic  area  to  fall  within  the  BCA.  This  could  mean  within  a  1000  km  radius  instead 
of  on  the  correct  one-third  of  the  world  for  an  Earth  Coverage  (EC)  antenna,  example. 

Spot  beams  also  allow  for  frequency  reuse,  which  allows  spectrum  to  be  separated 
geographically.  This  helps  prevent  unintentional  interference  from  neighboring  beams 
and  also  allows  the  finite  and  allocated  spectrum  to  be  used  several  times,  supporting 
more  users  (and  profit,  in  the  commercial  sector).  High  throughput  satellites  such  as 
INMARSAT  Global  Xpress  and  ViaSat  use  multiple  feeds  for  each  reflector  to  create 
patterns  of  spot  beams  to  allow  for  frequency  reuse.  03b  accomplishes  this  differently, 
taking  advantage  of  the  smaller  antenna  sizes  at  Ka-band  to  support  frequency  reuse  by 
allowing  for  many  steerable  antennas  to  be  placed  on  moderately  sized  satellites. 
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All  of  these  eharaeteristies  help  make  Ka-band  a  desired  frequeney  range  for  use, 
but  not  for  support  of  TDOA/FDOA  geoloeation.  TDOA/FDOA,  requiring  multiple 
satellites,  would  be  highly  unlikely  to  have  spot  beams  plaeed  in  the  same  loeation,  as 
one  of  the  benefits  is  that  they  ean  be  moved  to  support  dispersed  or  entirely  different 
users.  Narrow  beams  at  Ka-band  from  the  jammer  uplink  would  also  not  be  able  to 
simultaneously  target  multiple  satellites,  making  TDOA/FDOA  further  unlikely.  If 
multiple  beams  from  the  satellite  eould  be  steered  to  eover  the  same  area  to  perform 
the  geoloeation  that  would  remove  those  beams  from  their  primary  mission  of 
oommunieation  until  the  geoloeation  has  been  performed.  This  would  interrupt  serviee 
aeross  a  wider  area  than  was  originally  impaeted.  For  example,  every  time  a  jammer 
degraded  eommunieations  for  one  beam,  the  operator  would  then  interrupt 
eommunieations  on  other  beams  to  find  it.  If  an  adversary  realized  that  was  the  teehnique, 
it  eould  easily  eneourage  instead  of  diseourage  jamming.  This  also  presumes  that  the  user 
experieneing  interferenee  operates  the  other  beams  and  thus  moving  them  is  even  a 
ehoiee  (albeit  a  poor  one).  Leasing  a  seeond  or  third  beam  only  for  geoloeation  purposes 
is  far  from  ideal.  This  same  reasoning  is  applied  to  why  single  satellite/multiple  antenna 
geoloeation  methods  are  not  ideal  for  use,  partieularly  as  it  pertains  to  a  non-GEO 
oommunieations  eonstellation  (sueh  as  03b  for  the  foeus  of  this  thesis).  With  moving 
satellites  multiple  beams  would  need  to  be  leased  on  eaeh  satellite,  otherwise  the  ability 
to  perform  geoloeation  would  be  dependent  on  whieh  satellite  was  overhead  at  the 
moment.  Onee  again,  an  ideal  geoloeation  method  eould  use  a  single  spot  beam  on  a 
single  satellite  in  order  to  minimize  disruption  and  maximize  effieient  use  of  resourees. 

Single  antenna  methods  of  geoloeation  appear  ideal  for  use  by  03b  Networks 
to  support  interferenee  geoloeation  for  military  eustomers.  03b,  however,  does  not 
use  phased  array  antennas  that  are  eapable  of  measuring  phase  differenee  or  performing 
interferometry.  With  the  only  other  single  antenna  method  being  use  of  the  Doppler 
effeet,  it  appears  a  potentially  viable  method  for  use  that  requires  substantially 
fewer  resourees  than  other  eommon  methods.  This  will  be  addressed  in  greater  depth  in 
Chapter  III. 
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D.  CONCLUSION 


In  this  chapter,  it  was  established  that  geoloeation  of  EMI  emitters  is  an  important 
pieee  of  the  JSIR  proeess.  It  was  also  established  that  Ka-band  geoloeation  methods  must 
overcome  ehallenges  sueh  as  relatively  few  satellites  and  relatively  small  antenna 
footprints.  These  make  methods  of  geoloeation  sueh  as  TDOA  and  FDOA  unsuitable. 
Using  a  single  antenna  method  of  geoloeation  would  eonserve  assets,  and  use  of  Doppler 
shift  eould  be  a  way  to  effeetively  exeeute  it. 
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III.  03B  AND  SUITABILITY  FOR  GEOLOCATION 

TECHNIQUE(S) 


03b  is  a  satellite  eommunieations  company  that  provides  high-data-rate 
communications  to  customers  through  steerable  spot  beams.  03b  has  completed  various 
testing  objectives  with  the  U.S.  military,  among  those  a  test  with  7th  Fleet  during  Trident 
Warrior  2015,  work  with  the  Marine  Corps  Tactical  Systems  Support  Activity 
(MCTSSA),  and  demonstrations  with  the  Space  and  Naval  Warfare  Systems  Command 
(SPAWAR)  Common  Optical  Distribution  Architecture  (CODA)  lab.'^^ 

A.  CONSTELLATION  DESIGN 

The  03b  constellation  uses  a  circular  orbit  (<.001  eccentricity),  currently  in  one 
equatorial  (incline  of  about  .05°  or  less)  plane.  The  orbit  is  at  an  altitude  of  approximately 
8062  km,  which  is  a  relatively  low  (compared  to  Global  Positioning  System/GPS 
satellites)  MEO  orbit,  with  an  orbital  period  of  approximately  4.8  hours. The 
constellation  currently  contains  12  satellites. At  least  six  satellites  are  required  to 
provide  continuous  service  as  designed  (with  specified  overlap/make-before -break  hand- 
off). ^2  There  are  also  plans  to  expand  service  by  increasing  the  number  of  satellites  in  the 
current  plane,  and  there  has  been  some  hypothetical  discussion  about  eventually  adding 
two  inclined  planes  that  would  provide  polar  coverage. ^3 

B.  SATELLITE  CHARACTERISTICS 

03b  satellites  have  10  steerable  customer  antennas  and  two  steerable  gateway 
antennas. 34  The  antennas  are  maneuverable  to  provide  what  03b  calls  “optimal 

49  03b  Networks,  “03b  Government  USG  Bulletin,”  03bNetworks.com,  My  16,  2015, 
http://www.o3bnetworks.eom/wp-content/uploads/2015/08/USG-Bulletin_16JUL15.pdf. 

3®  Steven  H.  Blumenthal,  “Medium  Earth  Orbit  Ka  Band  Satellite  Communications  System,”  in 
MILCOM  2013  -  2013  IEEE  Military  Communications  Conference,  (San  Diego:  IEEE,  2013),  213-211. 

31  J.  J.  Shaw,  “03b  Networks  and  Why  Not?”  (see  chap.  1,  n.  20). 

32  Ibid. 

33  Ibid. 

34  Ibid. 
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coverage”  to  latitudes  up  to  45  degrees  North  or  South,  and  “extended  eoverage”  from 
45-62  degrees  North  or  South.^^  Eaeh  antenna  is  parabolic  with  a  circular  aperture,  the 
diameter  is  approximately  .4  m.^®  Satellites  currently  operate  in  commercial  Ka-band, 
which  uses  27-29  GHz  uplink  frequeneies  and  downlink  frequeneies  of  about  17-19 
GHz. 57  03b  generally  considers  the  widest  beamwidth  for  communieations,  “outer 
coverage,”  to  be  five  degrees  (2.5  degrees  off  boresight),  whieh  is  the  -6  dB  gain  point,  as 
can  be  seen  in  Figure  4.58  approximately  700  km  in  width  at  nadir,  as  shown  in 

Figure  5.59 


Figure  4.  03b  Satellite  Payload  Antenna  Reeeive  Pattern.®*’ 


—Rx  theta  cut 
— —  Rxphi  cut 


J.  J.  Shaw,  “03b  Networks  and  Why  Not?” 

5®  Ibid;  Ken  Mentasti,  email  message  to  author,  November  1 1,  2015. 

57  J.  J.  Shaw,  “03b  Networks  and  Why  Not?” 

58  Ibid. 

59  Ibid. 

®9  Source:  J.  J.  Shaw,  “03b  Networks  and  Why  Not?” 
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Figure  5.  03b  Beam  Coverage. 


C.  USEFUL  CHARACTERISTICS  FOR  GEOLOCATION 

As  was  covered  in  Chapter  II,  a  method  of  geolocation  of  EMI  sources  using 
a  single  antenna  would  be  ideal  as  it  would  minimize  the  impact  on  valuable 
communication  resources.  Without  a  phased-array  antenna  enabling  the  performance  of 
single-antenna  interferometry,  the  single  antenna  method  that  remains  is  use  of  the 
Doppler  effect,  although  the  exact  application  is  only  determined  through  the  rest  of  this 
thesis.  03b,  located  in  a  MEO  orbit,  has  a  greater  relative  velocity  to  a  stationary 
transmitter  as  compared  to  a  communications  satellite  in  a  GEO  orbit.  The  high 
frequencies  of  Ka-band  matched  with  the  greater  relative  velocities  also  contribute  to  a 
large  Doppler  shift — as  could  easily  be  seen  from  Equation  (1.2),  compared  to  other 
typical  frequency/orbit  combinations.  To  give  an  idea  of  the  magnitude  of  the  velocity 
difference,  an  03b  satellite  is  traveling  more  than  5.25  km/s,  whereas  a  GEO  satellite  is 
traveling  at  approximately  3.07  km/s.  That  is  a  difference  of  more  than  two  km/s,  or 


Source:  J.  J.  Shaw,  “03b  Networks  and  Why  Not?” 
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about  4877  miles  per  hour.  This  is  not  an  insignifieant  amount,  and  as  a  result  the 
Doppler  shift  is  mueh  larger,  easier  to  measure,  and  use. 

The  relative  motion  ereated  by  the  MEO  orbit  of  03b  and  the  ability  to  steer 
beams  is  also  important  for  use  with  a  method  that  would  employ  Doppler  shift 
measurement.  A  signifieant  reason  is  that  if  there  is  an  instant,  at  whieh  the  elosing 
veloeity  between  the  EMI  souree  and  satellite  is  zero,  then  the  Doppler  shift  would  also 
be  zero,  and  the  frequeney  reeeived  would  mateh  the  base  frequeney  of  the  EMI  souree. 
This  will  oeeur  when  the  interferenee  souree  passes  abeam  the  satellite’s  position  within 
the  beam  eoverage  area.  With  a  steerable  beam,  the  beam  ean  be  plaeed  along  a  line 
perpendieular  to  the  satellite’s  ground  traek  (abeam),  at  the  satellite’s  position  to  ensure 
this  oeeurs.  Eor  a  satellite  with  no  inelination  this  would  mean  the  same  meridian  as  the 
satellite  was  positioned  over.  This  is  also  made  possible  by  the  faet  that  the  satellite,  at 
MEO,  will  assuredly  pass  the  abeam  point  due  to  its  speed  at  faster  than  the  Earth’s 
rotation.  With  eommunieations  satellites  at  GEO,  motion  relative  to  the  Earth  is 
extremely  limited  by  eomparison  to  03b  and  in  many  eireumstanees  the  emitter  would 
never  pass  into  the  abeam  position.  This  is  a  key  faetor  in  using  Doppler  for  geoloeation, 
beeause  the  exaet  frequeney  of  an  interferenee  souree  eould  easily  be  unknown.  Without 
knowing  what  base  frequeney  is  added  to  the  Doppler  shift  (A/)  to  ereate  the  frequeney 
aetually  reeeived/measured  at  the  satellite,  there  is  no  way  to  determine  the  Doppler  shift 
present  at  the  varying  times  throughout  the  interferenee  event.  While  the  moment  the 
interferenee  source  passes  abeam  would  have  to  be  identified  somehow,  doing  so  would 
provide  an  LOP  perpendicular  to  the  satellite’s  motion  as  well  as  identify  the  base 
frequency.  Without  the  motion  due  to  03b ’s  unique  orbit  ensuring  the  emitter  passes 
abeam,  this  LOP  would  not  be  possible.  With  the  base  frequency  established  all  other 
received  frequency  data  could  now  be  used  for  the  Doppler  shift  it  contains  in  order  to 
correlate  a  position. 

The  steerable  beams  of  03b  satellites  also  enable  the  extended  coverage  area  of 
up  to  62  degrees  North  or  South,  potentially  providing  geoloeation  services  over  that 
entire  range.  A  constellation  of  12  satellites,  each  with  an  orbital  period  of  4.8  hours, 
further  means  that  a  satellite  is  passing  every  location  within  the  covered  latitudes  dozens 
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of  times  each  day  (four  times  per  satellite).  This  would  provide  multiple  opportunities  for 
an  EMI  source  to  be  located. 

D.  CONCLUSION 

The  03b  constellation  may  provide  a  useful  opportunity  for  a  commercially 
provided  geolocation  service.  03b ’s  satellites  and  constellation  differ  significantly  from 
most  communications  satellites  on  orbit  today  by  use  of  a  medium  Earth  orbit.  The 
velocity  of  the  satellites  relative  to  Earth  and  Ka-band  frequencies  serve  to  increase  the 
Doppler  shift  that  would  be  received.  The  number  of  satellites,  revisit  time,  and  steerable 
beams  all  work  to  increase  the  number  of  opportunities  and  the  responsiveness  of  any 
potential  geolocation.  These  characteristics  all  favorably  support  the  potential  ability  of 
03b  to  be  used  for  a  Doppler-based  geolocation  technique. 
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IV.  GEOLOCATION  TECHNIQUE 


Chapter  III  established  that  03b  has  unique  charaeteristies  that  make  possible  a 
new  method  of  geoloeation  using  the  reeeived  frequeney  from  an  emitter  over  time. 
While  the  eharaeteristies  of  Doppler  and  03b ’s  orbit  make  a  longitudinal  LOP  seem 
readily  aeeessible,  the  speeifies  of  that  teehnique  still  need  to  be  addressed,  as  does  a  way 
to  ereate  an  interseeting  LOP.  To  do  this,  this  ehapter  will  start  by  looking  at  another 
applieation  of  Doppler  shift,  and  then  adapt  it  to  a  simplified  geoloeation  problem.  Onee 
this  has  been  eompleted,  the  method  for  modeling  a  more  eomplex  situation  will  be 
introdueed.  This  model  was  then  used  throughout  the  final  seetion,  whieh  addresses 
added  eomplexity  and  realism. 

Geoloeation  of  an  EMI  souree  has  the  potential  to  be  a  eomplieated  problem  when 
attempting  to  eonsider  every  possible  situation  that  eould  oeeur  in  the  real  world.  These 
situations  eould  inelude  moving  emitters  at  speeds  ranging  from  those  of  aireraft  to  those 
of  boats  or  man,  various  altitudes  due  to  terrain  or  aireraft,  or  transmitters  that  vary  in 
frequeney,  among  many  more.  While  some  of  these  more  eomplieated  situations  have 
been  addressed  in  at  least  a  few  other  sourees,^^  this  thesis  was  not  intended  to  solve 
those  issues  but  rather  develop  and  provide  proof  of  eoneept  for  a  method  that  eould  be 
readily  applied  to  the  03b  Network  satellites  with  some  utility  for  DOD  operations.  To 
aeeomplish  that,  emitters  were  assumed  to  be;  stationary,  fixed  transmit  frequeney, 
eontinuously  emitting,  and  loeated  on  the  surfaee  of  a  rotating  Earth  as  represented  by  an 
ellipsoid.  Eurthermore,  the  03b  satellite  orbits  were  simplified  to  be  perfeetly  eireular 
(zero  elliptieity)  and  equatorial  (zero  inelination).  Also,  intentionally  left  out  is  the  small 
effeet  the  time  for  the  signal  to  travel  from  the  emitter  to  satellite  would  have  on 
aeouracy,  as  the  teehnique  laeks  the  preeision  to  make  that  important. 


Tina  L.  Chow,  “Passive  Emitter  Loeation  Using  Digital  Terrain  Data”  (master’s  thesis,  Binghamton 
University  of  State  University  of  New  York,  2001);  Hanna  Witzgall,  “A  Reliable  Doppler-Based  Solution 
for  Single  Sensor  Geoloeation,”  (paper  presented  at  the  IEEE  Aerospaee  Conferenee  ,  Big  Sky,  Mareh 
2013). 
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A.  TECHNIQUE  DEVELOPMENT 


Use  of  Doppler  shift  to  locate  or  refine  the  location  of  objects  by  aircraft/satellite 
is  not  new.  An  application  used  commonly  for  many  years  now  in  both  military  and 
civilian  applications  is  Doppler  beam  sharpening  of  radar  returns,  such  as  with  synthetic 
aperture  radar. ^3  All  that  is  necessary  then  to  use  Doppler  for  geolocation  of  EMI  sources 
instead  of  terrain  is  to  make  the  geolocation  situation  similar  to  Doppler  beam 
sharpening.  The  method  often  used  for  Doppler  beam  sharpening  contains  both  some 
differences  from  and  similarities  to  what  would  need  to  be  done  in  the  circumstances  of 
03b  locating  an  EMI  emitter.  In  the  use  of  Doppler  beam  sharpening  the  frequency  is 
known.  This  is  a  significant  advantage  over  locating  an  emitter  with  unknown  frequency 
as  the  Doppler  shift  can  be  calculated  directly  and  immediately  from  the  received 
frequency.  The  Doppler  shift  correlates  to  an  isodop,  which  is  a  hyperbola  shaped  EOF, 
anywhere  along  which  if  a  transmitter  (or  reflector)  was  positioned  would  cause  the  same 
Doppler  shift  to  occur.®^  Doppler  beam  sharpening  also  uses  a  pulse  emitted  by  the 
satellite  that  will  later  be  received  (by  the  satellite),  allowing  for  calculation  of  the  range 
from  the  amount  of  time  it  takes  the  signal  to  return.^^  is  also  a  significant  advantage 
as  it  provides  the  second  necessary  LOP  (a  circular  one)  necessary  to  limit  the  possible 
geolocation  to  four  points. Half  of  those  are  eliminated  by  the  known  factor  of  which 
side  the  radar  was  observing,  and  one  of  the  remaining  two  points  is  eliminated  by 
whether  the  Doppler  shift  was  positive  or  negative  (indicating  whether  the  emitter  lies  on 
the  symmetrical  isodops  either  ahead  or  behind  the  satellite’s  position). The  technique 
used  to  eliminate  positional  ambiguity  to  the  left  or  right  of  the  satellite’s  path  is  easily 
adapted  to  03b  as  the  beam  receiving  interference  can  be  steered  to  and  then  known  to  be 
located  on  one  side  or  the  other. 


^3  John  C.  Curlander  and  Robert  N.  McDonough,  Synthetic  Aperture  Radar  Systems  and  Signal 
Processing  (New  York:  John  Wiley  &  Sons,  1991),  17-26. 

64  Ibid.,  18-19. 

65  Ibid.,  18-19. 

66  Ibid.,  18-19. 

67  Ibid.,  18-19. 
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The  differences  highlighted  between  a  Doppler  beam  sharpening  application  and  a 
geolocation  of  EMI  sources  clarify  the  issues  that  need  to  be  resolved  in  order  to  use 
Doppler  shift  to  create  a  geolocation.  The  first  is  that  the  transmitted  frequency  needs  to 
be  identified  because  the  source  is  not  the  satellite,  as  it  is  in  Doppler  beam  sharpening. 
The  next  issue  is  how  to  obtain  a  range  LOP  from  Doppler  when,  once  again,  the  source 
of  the  transmission  is  not  the  satellite  and  therefore  the  origination  time  and  transit  time 
are  unknown.  Resolving  the  first  issue  would  yield  an  LOP  identifying  the  EMI  source  as 
abeam,  thereafter  requiring  an  LOP  giving  range  or  cross-track  distance  (the  second 
issue)  from  the  satellite  in  order  to  locate  the  EMI  source. 

To  evaluate  the  existing  issues,  it  is  useful  initially  to  simplify  the  problem  to  a 
flat,  non-rotating  Earth  and  observe  data  that  would  be  generated  and  usable  for  the 
geolocation  process,  and  how  that  data  appears  to  vary  with  location.  A  dataset  for  a  real 
world  situation  could  contain  the  received  frequency  (base  frequency  plus  Doppler  shift), 
which  is  easily  monitored,  and  the  timestamp  for  each  of  those  data  points.  A  satellite 
position  at  each  of  those  times  could  also  be  easily  obtained  through  use  of  GPS  in  the 
case  of  03b,  or  ephemeris  data  if  necessary,  both  of  which  can  have  very  good 
accuracy.  To  best  relate  the  data  received  to  a  transmitter  location,  MATLAB**’  was 
used  to  both  generate  and  plot  Doppler  shift  curves  with  varied  but  known  transmitter 
location  paths  (due  to  the  relative  motion  of  satellite/transmitter)  and  a  known  transmit 
(base)  frequency.  Once  the  relationship  between  the  Doppler  shift  over  time,  base 
frequency,  and  transmitter  location  over  time  are  characterized,  then  the  base  frequency 
and  location  can  become  variables  which  are  solved  for  by  analysis  of  the  received 
frequency  over  time. 

Consider  the  situation  shown  in  Eigure  6.  The  range,  R,  is  a  function  of  X,  Y,  and 
Z,  where  X  is  the  in-track  range,  Y  is  the  cross-track  range,  and  Z  is  the  satellite  altitude. 
The  in-track  range  varies  over  time  due  to  the  velocity  of  the  satellite  over  Earth,  Vsat,  and 
the  starting  in-track  range,  for  time  at  zero  seconds,  is  C.  The  dotted  red  lines  in  the 


Ken  Mentasti,  email  message  to  author,  November  1 1,  2015. 
MATLAB  is  a  registered  trademark  of  MathWorks. 
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figure  are  intended  to  show  the  change  over  time  from  the  initial  position.  The  range 
between  the  satellite  and  EMI  source  at  any  point  in  time  is  given  by  Equation  (1.5). 


=  +  +  (1-5) 

From  this  equation,  the  first  derivative  can  then  be  easily  found  to  give  the 
relative  velocity  necessary  for  the  Doppler  shift  equation,  Equation  (1.2).  When  added  to 
the  base  frequency,  this  results  in  the  frequency  that  would  actually  be  received  by  a 
satellite. 


Figure  6.  Flat  Non-rotating  Earth  Geolocation  Geometry. 


c 


Figure  7  was  created  in  MATEAB.  The  figure  shows  the  frequency  that  would  be 
received  (Y-axis  of  Figure  7)  by  an  03b  satellite  traveling  at  orbital  velocity  over  a  fiat, 
non-rotating  Earth,  for  various  in-track  distances  (X-axis,  positive  is  ahead  of  satellite) 
and  cross-track  distances  (each  of  four  lines  plotted,  better  resolution  in  Figures  7-9)  of 
an  emitter  with  a  base  frequency  of  30  GHz.  When  related  to  Figure  6,  Figure  7 
represents  a  situation  where  in-track  distance  is  based  on  Vsat=5.255  km/s  and  C=1000 
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km.  Z  would  be  an  altitude  of  8062  km  and  each  curve  in  Figure  7  would  represent  a 
different  Y,  700  km,  800  km,  900  km,  and  1000  km. 


Figure  7.  Received  Frequency  as  a  Function  of  Emitter  Relative  Position. 


While  Figure  7  gives  some  perspective,  it  is  not  precise  enough  to  be  useful.  For 
that,  closer  observation  is  needed  at  three  general  areas  of  concern:  prior-to-abeam, 
abeam,  post-abeam.  Figure  8  is  a  higher  resolution  graph  for  when  the  four  emitters  are 
approximately  1000  km  ahead  of  the  receiving  satellite  (and  closing)  at  their  respective 
cross  track  distances  (given  in  the  legend).  Figure  9  is  a  higher  resolution  graph  as  the 
emitters  transit  the  abeam  location  at  their  respective  cross  track  distances,  from  .5  km 
ahead  of  to  .5  km  behind  the  satellite.  In  Figure  9  the  curves  intersect  at  x=0  (abeam)  and 
are  indistinguishable  from  each  other  at  the  scale  given.  Figure  10  portrays  the  received 
frequency  as  the  emitters  are  approximately  1000  km  behind  the  satellite  (and  getting 
farther). 
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Figure  8.  Doppler  Shift  from  Emitters  Located  ahead  of  Satellite. 


In  Track  Distance  (km) 


Figure  9.  Doppler  Shift  from  Emitters  Transiting  abeam  Satellite. 


In  Track  Distance  (km) 
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Figure  10.  Doppler  Shift  from  Emitters  Located  behind  Satellite. 


In  Figures  8-10  several  useful  pieces  of  information  can  be  readily  observed  that 

demonstrate  analysis  of  the  Doppler  shift  Equation  (1.2).  The  first  is  that,  from  Figure  8 

and  Figure  10,  at  locations  where  the  emitter  is  not  abeam,  there  is  separation  in  the 

Doppler  shifts  that  can  be  related  to  the  cross-track  distance.  This  readily  correlates  with 

the  Doppler  equation,  as  relative  velocity  between  the  transmitter  and  emitter  is  greater 

when  the  satellite  velocity  aligns  closer  with  the  relative  position  vector.  This  means  that 

given  emitters  at  varying  cross-track  distances  and  equal  in-track  distances,  the  closest 

emitter  has  the  greatest  magnitude  of  Doppler  shift  at  any  in-track  position  except  at 

abeam.  This  can  be  identified  in  Figure  8  by  the  700  km  range  curve  having  more 

positive  Doppler  shift  than  the  other  ranges  graphed,  and  in  Figure  10  by  the  700  km 

range  curve  having  more  negative  Doppler  shift  than  the  other  ranges  displayed.  At 

abeam,  there  is  no  component  of  the  satellite  velocity  vector  aligned  with  the  relative 

position  vector  as  they  are  perpendicular,  causing  there  to  be  no  Doppler  shift  regardless 

of  cross-track  distance.  The  rate  at  which  the  relative  velocity  changes  is  also  greatest  at 

the  abeam  point.  This  can  be  evidenced  by  the  steeper  slope  of  the  Doppler  shift  curve  in 

Figure  9,  which  displays  the  data  in  close  proximity  to  the  abeam  point  as  opposed  to  the 

lesser  magnitude  slopes  at  locations  farther  away  from  the  abeam  point  in  Figures  8  and 
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10.  Figures  8-10  all  display  a  one  km  in-traek  seetion  in  order  to  make  this  eomparison. 
From  this  information,  the  eonelusion  ean  be  drawn  that  the  maximum  slope  and  an 
infleetion  point  oeeur  at  the  abeam  loeation. 

While  the  eoneepts  ean  be  demonstrated  impreeisely  from  the  figures,  more 
preeise  information  ean  be  obtained  by  using  a  little  math.  When  the  data  are  evaluated 
they  will  be  in  the  form  of  points,  taken  at  a  given  sample  rate.  The  first  step  to  make  the 
data  easier  to  use  would  be  to  fit  a  eurve  to  the  data.  This  is  easily  aeeomplished  through 
use  of  any  number  of  programs;  MATLAB  and  a  third  order  eurve  fit  were  used  for  this 
thesis.  The  next  pieee  of  information  that  is  important  to  obtain  is  the  loeation  of  the 
maximum  slope,  or  the  infleetion  point,  whieh  are  eollocated  in  the  graph.  This  ean  be 
done  by  finding  the  root  of  the  seeond  derivative  of  the  frequeney  eurve  (where  the  first 
derivative  is  the  slope).  The  solution  of  this  is  the  time  at  whieh  the  emitter  was  abeam 
the  satellite.  Assuming  that  the  satellite  is  at  a  known  loeation  with  whieh  the  time  ean  be 
assoeiated,  the  first  of  two  neeessary  LOPs  for  the  geoloeation  ean  be  established.  The 
frequeney  reeeived  at  the  abeam  time  is  also  known  to  be  the  base  frequeney  sinee  there 
is  no  Doppler  shift  at  that  loeation,  as  previously  diseussed. 

The  seeond  LOP  will  establish  the  eross-traek  distanee.  Now  that  the  base 
frequeney  has  been  ealeulated,  eurves  ean  be  ereated  whieh  illustrate  potential  reeeived 
frequeney  data  from  simulated  emitters  with  identieal  base  frequeney  to  the  jammer  and 
loeated  at  varying  eross-traek  distanees,  similar  to  what  was  done  in  Figure  7.  The 
differenee  is  that  the  eurves  ean  be  tailored  to  the  same  extent  of  time  that  the  emitter 
frequency  data  has  been  collected  for  and  plotted  simultaneously  against  the  collected 
data,  with  the  curves  positioned  so  that  the  base  frequency  value  on  the  curve  will  occur 
at  the  abeam  time  (the  curves  will  all  intersect  here  because  of  this,  similar  to  Figure  9). 
For  the  same  scenario  as  in  Figures  8-10,  a  received  frequency  curve  for  a  jammer  can  be 
simulated  at  an  additional  cross  track  distance  (840  km  is  used  for  the  jammer  in  Figure 
11).  The  received  frequency  plots  correlating  to  when  the  emitter  was  approximately 
1000  km  in  front  of  the  satellite  are  shown  in  Figure  1 1 . 
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Figure  11.  Received  Frequency  from  EMI  and  Simulated  Emitters. 


In-Track  Distance  (km) 

While  the  point  at  which  the  cross-track  distance  will  be  interpolated  from  for  the 
second  EOF  does  not  need  to  be  at  the  farthest  distance/time  from  the  abeam  point,  as  in 
Figure  11,  that  is  where  the  most  separation  in  frequency  occurs  between  the  curves.  The 
cross-track  distances  cannot  be  compared  at  the  abeam  point  as  all  the  curves  intersect 
there,  and  with  no  separation  in  the  frequencies  the  interpolation  cannot  be  performed. 
Interpolation  of  a  cross-track  distance  for  the  jammer  is  easily  done  and  the  results  are 
almost  obvious  from  Figure  1 1  already.  For  a  fixed  in-track  distance  (or  time)  data  points 
correlating  to  the  simulated  emitters  (at  known  cross-track  distances)  from  the  graph  are 
pulled  that  consist  of  the  frequency  received  on  that  curve  at  the  comparison 
distance/time  as  the  x-value  and  the  cross-track  distance  of  that  curve  (given  in  the 
legend)  as  a  y-value.  These  values  are  then  plotted  and  a  curve  fitted  to  them,  which  is 
Figure  12.  From  this  point  the  received  frequency  for  the  jammer  at  the  comparison 
distance/time  becomes  the  x-value,  and  the  equation  of  the  best  fit  line  will  return  the 
cross  track  distance. 
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Figure  12.  Graph  for  Interpolation  of  Cross-Track  LOP. 
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While  the  simplified  execution  of  the  method  for  a  flat,  non-rotating  Earth  and  use 
of  a  simulated  jammer  can  be  shown  to  work,  a  more  realistic  situation  creates  numerous 
challenges  that  need  to  be  resolved. 


B.  MODELING 


Before  addressing  adding  complexity  to  the  technique,  it  was  necessary  to 
develop  a  model  that  would  enable  accuracy  to  be  checked  so  that  progress  in  dealing 
with  the  added  complexity  could  be  tracked.  This  model  would  ensure  that  data  generated 
for  the  EMI  emitter  would  accurately  reflect  real-world  data,  and  also  that  solutions  were 
not  inadvertently  input  to  the  MATLAB  program  that  was  created  to  perform  the 
geolocation.  The  modeling  was  done  through  the  use  of  Systems  Tool  Kit  (STK),  a 
product  of  Analytical  Graphics,  Inc.,  which  allows  for  highly  realistic  modeling  and 
simulations  to  be  created  and  analyzed  for  a  wide  variety  of  situations.  STK  is  used  as  an 
analysis  tool  by  both  the  civil  and  military  space  communities  and  is  accepted  as  an 
industry  standard.  In  this  instance,  it  was  used  to  analyze  the  communications  link 
between  Earth  and  space,  specifically  emitters/jammers  and  a  modeled  03b  satellite. 
Ground-based  emitters  were  placed  at  various  locations  to  represent  the  majority  of  the 
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03b  coverage  area  (when  considering  results  for  north  and  south  latitudes  are 
symmetrieal)  in  order  to  provide  some  insight  to  performanee  over  varying 
latitudes/ranges.  An  03b  satellite  model  was  also  ereated,  with  the  orbit  simplified  for 
ease  of  use  as  previously  diseussed.  The  frequency  used  for  the  link  was  30  GHz,  which 
is  slightly  higher  than  03b’s  commereial  Ka-band,  and  falls  within  the  military  Ka-band. 

STK  is  able  to  generate  a  detailed  report  on  the  link  budget  over  time.  For 
modeling  the  situation  for  this  thesis,  the  data  used  was  restrieted  to  what  would  be 
observable  by  the  satellite  if  the  loeation  of  the  emitter  was  not  known  (for  example,  no 
range  data  for  the  link  was  used).  The  data  used  consisted  of  the  frequeney  received  by 
the  satellite  from  the  emitter,  and  a  time  stamp  for  each  frequency;  this  was  then  exported 
to  Mierosoft  Excel.  In  Microsoft  Excel,  the  time  of  each  received  frequency  was 
modified  to  be  a  simple  “elapsed  time”  that  simply  started  at  the  first  point  with  zero  and 
eounted  upward  at  an  interval  matching  the  sample  rate  (which  would  be  known,  since 
that  is  a  funetion  of  the  satellite  monitoring  and  not  the  emitter).  The  MATLAB  program 
would  then  import  the  STK-generated  data  from  Microsoft  Excel  for 
proeessing/geoloeation.  Also  entered  into  the  program  was  the  satellite’s  longitude  at  the 
start  of  the  STK  scenario.  After  the  program  computed  the  geolocation  estimate  for  the 
emitter,  this  was  then  compared  to  the  actual  location  of  the  emitter  in  STK  to  evaluate 
accuracy. 

C.  ADDING  COMPLEXITY 

Now  that  a  method  and  model  have  been  established  to  track  how  well 
eomplexity  is  addressed,  it  is  time  to  increase  the  realism  from  the  fiat,  non-rotating 
Earth,  to  the  real  world,  which  is  not  fiat  (or  spherical),  and  is  rotating.  While  the  latitude 
portion  of  the  geolocation  solution  will  change  signifieantly  when  adding  complexity,  the 
other  EOF  (longitude),  must  be  eomputed  first  and  ehanges  little  from  the  previously 
diseussed  flat,  non-rotating  Earth  scenario.  The  data  used  for  the  geoloeation  will  consist 
of  a  time  and  a  reeeived  frequeney,  like  the  data  shown  in  Table  1.  A  eomplete  set  of  data 
for  one  modeled  emitter  ean  be  found  in  Appendix  A.  The  middle  column  of  Table  1 
simply  counts  the  length  of  time  interferenee  has  been  reeeived,  which  is  an  easier  format 
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of  time  to  manipulate/use  in  the  MATLAB  program  than  UTCG.  With  a  non-flat  Earth, 
distanee  and  time  are  no  longer  related  as  simply  as  in  the  previous  seetion,  beeause 
translational  veloeity  of  the  emitter  varies  with  latitude,  whieh  is  unknown.  This  makes  it 
neeessary  to  work  with  the  data  as  it  is  refereneed  to  time. 


Table  1 .  Sample  Data  Exeerpt  for  Geoloeation. 


Time  (UTCG) 

Recv  Interference  Elapsed  Time 

Rcvd.  Frequency  (Hz) 

29  Mar  2016  20:11:46.000 

12 

30000042969.19020 

29  Mar  2016  20:11:47.000 

13 

30000042884.06240 

29  Mar  2016  20:11:48.000 

14 

30000042798.91810 

29  Mar  2016  20:11:49.000 

15 

30000042713.75750 

29  Mar  2016  20:11:50.000 

16 

30000042628.58040 

29  Mar  2016  20:11:51.000 

17 

30000042543.38700 

29  Mar  2016  20:11:52.000 

18 

30000042458.17730 

The  reeeived  frequency  data  are  fitted  with  a  curve  which  is  then  plotted,  much 
like  Figure  9,  with  the  frequency  on  the  y-axis,  but  now  with  time  on  the  x-axis  instead  of 
distance,  as  shown  in  Figure  13.  The  inflection  point  of  the  curve  is  then  found  from  the 
second  derivative,  as  was  done  in  the  previous  scenario.  The  time  at  which  this  inflection 
point  occurs  is  once  again  the  abeam  point,  and  the  received  frequency  at  that  time  is  the 
base  frequency,  all  based  on  the  previously  established  concepts.  The  base  frequency  will 
be  used  later  for  the  comparison  curve  generation  to  find  the  latitude  EOF  (same  process 
as  the  cross-track  EOF  previously).  With  the  abeam  time  for  the  emitter  identified,  the 
longitude  of  the  satellite  at  that  time  (known  to  the  satellite  operator  or  user)  is  equal  to 
the  longitude  of  the  emitter,  and  the  first  EOF  is  complete.  The  times  for  the  data  are  then 
adjusted  so  that  the  abeam  time  is  set  to  zero  seconds,  and  positive  time  represents  time 
prior  to  abeam,  while  negative  time  will  represent  data  for  when  the  emitter  is  past  the 
abeam  point.  This  adjustment  is  somewhat  arbitrary,  but  simplifies  comparison  to  the 
generated  curves  for  the  latitude  EOF  process  as  the  y-intercept  for  all  curves  will  be  the 
base  frequency.  After  this  adjustment  has  been  completed  the  resultant  graph  is  shown  in 
Figure  14. 
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Figure  13.  Received  Frequency  Plot  for  EMI  Emitter  from  STK  Data. 


Figure  14.  Received  Frequency  Plot  Adjusted  after  Abeam  Time  Identified. 
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The  next  portion  of  the  geolocation  is  the  calculation  of  the  latitude.  When  the 
flat,  non-rotating  Earth  is  actually  an  elliptical,  rotating  Earth,  and  theoretical  emitter 
location  could  drastically  vary  in  latitude,  the  situation  is  significantly  more  complex  than 

with  a  flat  Earth.  It  is,  however,  necessary  complexity  in  order  to  obtain  useful  results. 
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Received  Freq  (Hz)  vs.  Time 


The  accuracy  of  the  latitude  solution  depends  on  the  ability  to  accurately  simulate  the 
received  frequency  curve  of  simulated  emitters  from  various,  but  known  locations.  This 
will  then  be  compared  to  data  from  the  emitter  of  interest,  specifically  a  graph  similar  to 
Figure  14.  This  will  only  produce  beneficial  results  if  the  internally  generated  data  for 
comparison  can  mimic  what  would  actually  be  produced  if  an  emitter  was  at  that  location 
on  the  Earth,  which  means  that  the  Earth’s  shape  and  rotation  must  be  accounted  for. 

The  first  issue  is  that  a  different  reference  system  is  needed  since  the  Earth’s 
surface  is  no  longer  considered  fiat.  Additionally  complicating  this  issue  is  that  the  Earth 
is  also  not  a  sphere  (the  polar  axis  is  shorter  than  an  equatorial  axis),  and  considering  it  as 
such  would  negatively  affect  results,  particularly  closer  to  the  poles  and  locations  where 
the  radius  would  not  match  that  of  a  spherically  modeled  Earth.  The  next  issue  is  that  the 
Earth  is  rotating,  and  the  rotating  Earth  results  in  different  velocities  of  stationary 
(relative  to  Earth)  emitters  at  differing  latitudes.  This  is  due  to  the  different  circumference 
lengths  (when  the  radii  of  the  circle  is  measured  perpendicular  to  the  polar  axis)  of 
different  latitudes  while  all  points  on  the  Earth  have  the  same  angular  velocity.  The 
velocity  of  emitters  varying  with  latitude  matters  because  Doppler  shift  is  directly  related 
to  the  relative  velocity  between  the  satellite  and  emitter. 

Eor  ease  of  use,  a  Cartesian  reference  system  is  desirable.  In  this  situation,  it  will 
be  based  at  the  center  of  the  Earth — because  that  point  does  not  depend  on  the  Earth’s 
shape.  Eor  simplification,  the  satellite  will  be  considered  fixed  in  the  reference  system 
and  only  the  emitter  will  move,  representing  all  relative  motion.  This  reference  system 
will  consist  of  X  in  the  equatorial  plane  through  the  center  of  the  Earth  to  the  satellite,  Y 
perpendicular  to  X  in  the  equatorial  plane,  and  Z  will  coincide  with  the  polar  axis  of  the 
Earth  reference  ellipsoid.  Eor  the  satellite,  which  operates  at  a  fixed  altitude  (and  thus 
semi-major  axis,  which  will  be  the  “x”  coordinate),  and  is  at  a  constant  zero  “z” 
coordinate  due  to  the  equatorial  orbit,  the  Cartesian  coordinates  are  simple.  Using  this 
reference  system  also  means  the  “y”  coordinate  is  zero,  and  Earth  referenced  longitude 
does  not  matter — as  long  as  the  satellite’s  Earth  referenced  latitude  is  known  at  some 
point  during  the  data  collection  the  change  in  longitude  from  the  known  point  can  be 
tracked,  for  purposes  of  identifying  the  emitter’s  abeam  position  later.  It  is  also  worth 
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noting  that  these  coordinates  are  constant  over  time.  This  reference  system  is  illustrated 
in  Figure  15. 


Figure  15.  Reference  System. 


With  the  satellite’s  Cartesian  position  established,  the  Cartesian  position  of 
simulated  emitters  must  now  be  calculated  in  order  to  determine  their  relative  positions. 
The  emitter’s  relative  position  in  this  reference  system  changes  with  time,  and  because  of 
this  velocity  becomes  a  consideration.  The  emitter  is  also  on  the  surface  of  the  ellipsoidal 
Earth  and  received  frequency  curves  must  be  calculated  for  varying  latitudes, 
complicating  the  computation  of  coordinates  for  the  emitter  more  than  the  satellite.  To 
help  address  the  ellipsoidal  Earth  issue,  instead  of  creating  comparison  curves  for 
simulated  emitters,  similar  to  Eigure  11,  identified  by  cross-track  distance  in  km  from  the 
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nadir  point  of  the  satellite  (700  km,  800  km,  ete.),  the  eomparison  eurves  will  be 
refereneed  to  geodetie  latitude.  In  order  to  use  geodetie  latitude  as  a  referenee  while 
performing  the  math  to  generate  the  data  with  a  Cartesian  eoordinate  system,  eonversion 
is  neeessary.  Fortunately,  this  proeedure  is  detailed  in  SMADP^  The  first  step  is  to 
eompute  the  “ellipsoidal  radius  of  eurvature  in  the  meridian,”  shown  in  Equation  (1.5), 
using  the  semi-major  axis  of  the  ellipsoid  (a)  and  the  eeeentrieity  of  the  ellipsoid,  {e)J^ 
This  is  then  used  to  eomplete  the  eonversion  in  Equation  (1.5).  The  additional  inputs  are 
“ellipsoidal  height”  {h),  whieh  for  an  emitter  on  the  surfaee  of  Earth  will  be  zero,  the 
geodetie  latitude  {</)),  which  is  known  for  the  curves  being  generated,  and  longitude  (/I), 
which  will  be  explained  shortly.  ^2 

=  a  I  ^/l--e^'sin^^  [see  footnote]  23  (1-5) 


i)cos^/icos/l 

Y 

= 

K+' 

It]  cos  sin /I 

[see  footnote]  24 

(1.5) 

Z 

((1-" 

)A^+/t)sin^ 

The  most  useful  longitude  for  this  procedure  is  relative  to  the  satellite’s  longitude, 
as  the  goal  at  this  point  is  to  generate  simulated  received  frequency  data  which  includes 
Doppler  shift  that  varies  with  relative  position.  Eongitude  of  the  emitter  relative  to  the 
satellite  changes  with  time  as  a  result  of  both  the  rotation  of  the  Earth  and  the  orbit  of  the 
satellite.  A  useful  way  to  relate  the  longitude  to  time  is  to  use  rotational  velocity,  which  is 
not  impacted  by  the  latitude  of  the  emitter  as  translational  velocity  would  be.  The 
rotational  velocity  of  the  satellite  relative  to  the  emitter  is  the  difference  of  the  rotational 
velocity  of  Earth  and  the  rotational  velocity  of  the  satellite  in  orbit.  The  difference  in 
rotation  velocity  coupled  with  time  provides  the  required  relative  longitude  equation.  It  is 


20  Larson  and  Wertz,  Space  Mission  Analysis  and  Design,  900  (see  chap.  1,  n.  7). 

21  Ibid. 

22  Ibid. 

23  Ibid. 

24  Ibid. 
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important  that  the  equations  are  variable  with  time  and  no  other  parameters,  as  time  is 
one  of  only  two  data  types  that  result  from  and  ean  be  direetly  eompared  to  the  actual 
emitter  (the  other  being  frequency).  The  way  this  technique  is  applied  is  that  the  relative 
rotational  velocity  {vrot)  in  radians/second  is  calculated  as  shown  in  Equation  (1.5),  where 
4  is  the  satellite’s  orbital  period  and  4  is  Earth’s  rotational  period  (the  length  of  a  sidereal 
day). 


V  = 

rot 


f— 1 

\  J 

(1.5) 


Equation  (1.5)  does  not  vary  with  time,  and  becomes  an  input  for  solving  the 
relative  longitude  as  a  function  of  time.  This  is  done  by  multiplying  Vrot  by  the  amount  of 
time  ahead  or  behind  the  abeam  position  (which  has  already  been  calculated)  to  output 
the  relative  longitude  (/I)  for  conversion  to  Cartesian  coordinates,  as  shown  in  Equation 
(1.5).  A  quick  check  for  this  is  that  at  time  {t)  =  0  seconds,  the  relative  longitude  would 
also  be  zero,  and  therefore  be  on  the  same  meridian  as  the  satellite. 

^  =  v^J  (1.5) 


With  both  the  satellite  and  the  simulated  emitter’s  locations  now  known  as  a 
function  of  time,  the  relative  velocity  between  the  two  can  be  calculated.  Eor  this,  the 
range  will  first  be  calculated  using  Equation  (1.5).  This  is  also  shown  in  Eigure  15, 
denoted  by  “i?.”  In  the  program  created,  the  change  in  range  between  the  satellite  and 
emitter  over  a  small  increment  of  time  (ranges  Ri  and  R2)  was  used  to  obtain  the  relative 
velocity,  as  shown  in  Equation  (1.5).  This  velocity  then  used  to  obtain  the  Doppler  shifts 
and  added  to  the  base  frequency  (/),  which  was  found  at  the  abeam  position,  to  obtain 
received  frequency  (free)  curves  for  the  simulated  emitters,  as  shown  in  Equation  (1.5). 

R  =  +(T.  -yef  +(A  (1-5) 


rel 


Mi  _{R,-Rf) 
At  At 


(1.5) 


free  -  f  +  ^f  -  f  +—'^rel 

c 


(1.5) 
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In  this  instance,  received  frequency  curves  were  generated  for  the  simulated 
emitters  at  every  .5  degrees  of  latitude  from  3.5  degrees  below  to  3.5  degrees  above  beam 
center  (the  latitude  of  the  beam  center  would  be  known  by  the  operator).  This 
approximates  latitudes  that  would  be  within  the  BCA,  in  order  to  lessen  the  programming 
burden.  This  does  not  take  into  account  changes  in  BCA  with  latitude,  but  that  was  not 
necessary  to  prove  the  concept  as  long  as  the  method  could  locate  emitters  away  from 
beam  center.  Repeating  the  process,  curves  could  be  calculated  within  the  BCA  at  as 
small  an  interval  as  desired.  Adding  curves  to  represent  locations  farther  from  beam 
center  would  accommodate  a  larger  BCA.  For  the  geolocation  to  be  effective,  the 
comparison  curves  just  need  to  encompass  the  possible  latitudes  of  the  emitter,  as  those 
create  the  points  a  curve  will  be  fit  to  for  the  interpolation.  The  comparison  curves 
generated  for  this  instance  were  based  on  a  beam  centered  at  57.75°  N  and  are  shown  in 
Figure  16.  The  conditions  demonstrated  could  represent  the  Foresight  of  the  beam  aimed 
at  a  user  located  at  57.75°  N  and  the  emitter  of  interest  somewhat  further  from  the 
Foresight  but  still  within  the  BCA.  This  was  done  to  demonstrate  that  the  emitter  of 
interest  (located  at  60°  N)  does  not  need  to  be  (and  presumably  often  would  not  be) 
located  perfectly  at  the  center  of  the  BCA  for  the  geolocation  method  to  still  be  effective. 

The  next  step  is  to  compare  the  data  from  the  emitter  of  interest  (jammer)  in 
Figure  14  to  the  comparison  curves,  as  was  discussed  and  done  for  Figure  11  in  the 
previous  section.  Figure  17  shows  the  plot  of  these  curves  together,  focused  on  the  end  of 
the  curves  where  the  separation  is  most  apparent.  From  the  same  time  value  for  each  of 
the  simulated  emitter  curves  (at  known  latitudes),  data  points  are  taken  that  consist  of  the 
frequency  received  on  that  curve  at  the  comparison  time  as  the  x-value,  and  the  latitude 
of  the  emitter  for  that  curve  as  a  y-value.  Figure  18  displays  the  curve  created  from  these 
points  (latitude  is  in  radians  as  the  program  created  uses  radians  for  most  intermediate 
steps).  The  received  frequency  for  the  EMI  emitter  at  the  comparison  time  is  input  as  an 
x-value,  and  the  equation  of  the  best  fit  line  returns  the  latitude.  Results  will  be  discussed 
in  the  following  chapter. 
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Received  Frequency  (Hz) 


Figure  16.  Received  Frequency  Curves  of  Simulated  Emitters  at  Various 

Latitude. 


Figure  17.  Comparison  of  Received  Frequency  between  EMI  and  Simulated 

Emitters. 
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Figure  18.  Graph  for  Interpolation  of  Latitude  of  EMI  Emitter. 


D.  CONCLUSION 

This  chapter  developed  a  technique  to  take  advantage  of  data  that  could  be 
obtained  and  used  due  to  the  unique  characteristics  of  03b’s  satellites.  The  result  of  this 
technique  is  the  MATEAB  program  in  Appendix  B  which  is  capable  of  outputting  a 
geolocation  estimate  for  input  data  generated  by  STK,  which  models  realistic  conditions 
with  consideration  to  the  constraints  and  assumptions  addressed  at  the  beginning  of  the 
chapter.  The  next  chapter  will  address  the  results  obtained  from  this  geolocation 
technique  by  comparison  of  the  estimated  locations  of  emitters  to  the  actual  locations. 
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V.  RESULTS,  FUTURE  WORK,  AND  CONCLUSIONS 


Now  that  the  technique  for  the  new  geolocation  method  has  been  developed,  the 
results  must  be  evaluated  for  accuracy  and  then  analyzed  to  determine  what  application 
or  utility  they  serve  in  the  current  form.  For  evaluation  of  accuracy  with  this  method 
latitude  and  longitude  may  be  analyzed  independently  because  different  procedures  were 
used  to  create  their  respective  LOPs.  While  a  detailed  error  analysis  was  not  conducted, 
some  possible  sources  of  error  are  identified.  Future  work  to  build  on  this  thesis  will  also 
be  discussed  before  concluding  this  chapter  and  this  thesis. 

A.  RESULTS 

In  order  to  characterize  the  accuracy  of  this  method  through  a  wide  range  of 
possible  locations,  emitters  at  15°N  000°E,  30°N  000°E,  45°N  000°E  and  60°N  000°E 
were  modeled  in  STK  and  then  a  geolocation  estimate  was  computed  for  each  emitter 
using  a  MATLAB  program  (see  Appendix  B).  Eor  each  of  the  four  emitters  an  estimated 
location  was  computed  in  three  circumstances — one  with  the  emitter  at  the  beam  center, 
one  with  the  beam  centered  at  three  degrees  latitude  below  the  emitter,  and  again  with  the 
beam  centered  at  three  degrees  latitude  above  the  emitter.  Three  degrees  latitude  from 
Foresight  was  chosen  because  the  emitter  must  be  within  the  BCA,  and  the  MATLAB 
program  computes  comparison  curves  out  to  3.5  degrees  off  Foresight  to  simulate  the 
BCA.  Three  degrees  off  Foresight  therefore  represents  a  situation  where  the  emitter 
would  be  toward  the  edge  of  the  BCA,  but  still  within  it.  To  simulate  the  knowledge  a 
satellite  operator  would  have  of  the  satellite  beam’s  Foresight  location,  and  to  test 
accuracy  for  emitter  locations  within  the  BCA  but  not  at  the  Foresight  location,  the 
program  includes  an  input  for  centering  the  beam  at  specified  latitude.  This  can  be  found 
in  Section  1  of  the  MATLAB  program  (Appendix  B).  Section  1  of  the  program  also 
includes  an  input  for  the  actual  emitter  location,  which  would  obviously  not  be  known 
under  normal  circumstances,  but  was  used  in  order  to  calculate  error  distances  of  the 
estimated  emitter  location  from  the  actual  emitter  location,  performed  in  Section  8  of  the 
MATLAB  program. 
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Tables  2-4  display  the  estimated  emitter  loeations  ealeulated  using  the 
geoloeation  teehnique,  and  the  distanee  (km)  in  error  that  the  estimated  positions  were 
from  the  aetual  emitter  loeations  for  a  situation  where  the  emitter  falls  at  the  beam’s 
boresight  as  well  as  where  it  would  be  at  the  northern  or  southern  edges  of  the  BCA.  The 
distanee  used  for  the  error  in  these  tables  was  ealeulated  by  MATLAB  for  the  differenee 
in  eoordinates  between  the  positions,  over  the  surfaee  of  the  Earth  as  represented  by  an 
ellipsoid. 


Table  2.  Aetual  and  Estimated  Emitter  Eoeations — Emitter  Eoeated  at  Beam’s 

Boresight  Eatitude. 


Aetual  Emitter  Position 
(DD.DDDD) 

Estimated  Emitter  Position 
(DD.DDDD) 

Error  Distanee  (km) 

15°N  000°E 

15.0513°N  000.0880°W 

11.040 

30°N  000°E 

29.9960°N  000.1175°W 

11.348 

45°N  000°E 

45.0141°N  000.1616°W 

12.837 

60°N  000°E 

60.0445°N  000.2455°W 

14.557 

Table  3.  Aetual  and  Estimated  Emitter  Eoeations — Beam’s  Boresight  3°S  of 

Emitter  Eatitude. 


Aetual  Emitter  Position 
(DD.DDDD) 

Estimated  Emitter  Position 
(DD.DDDD) 

Error  Distanee  (km) 

15°N  000°E 

14.9394°N  000.0880°W 

11.601 

30°N  000°E 

29.9838°N  000.1175°W 

11.481 

45°N  000°E 

45.0246°N  000.1616°W 

13.030 

60°N  000°E 

60.0445°N  000.2455°W 

14.557 
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Table  4.  Actual  and  Estimated  Emitter  Eocations — Beam’s  Boresight  3°N  of 

Emitter  Eatitude. 


Actual  Emitter  Position 
(DD.DDDD) 

Estimated  Emitter  Position 
(DD.DDDD) 

Error  Distance  (km) 

15°N  000°E 

14.9394°N  000.0880°W 

11.601 

30°N  000°E 

30.0082°N  000.1175°W 

11.375 

45°N  000°E 

45.0211°N  000.1616°W 

12.955 

60°N  000°E 

60.0585°N  000.2455°W 

15.156 

These  results,  for  latitudes  from  15°N  to  60°N  (it  would  be  symmetrical  for  the 
South)  and  varying  emitter  positions  within  the  BCA  give  a  maximum  error  distance  of 
15.505  km.  One  noticeable  piece  of  information  from  these  three  tables  is  that  the 
longitude  error  for  each  emitter,  at  each  of  the  three  relative  emitter  positions  within  the 
BCA  is  the  same.  Eor  example,  the  estimated  longitude  for  the  emitter  actually  located  at 
15°N  000°E,  is  000.0880°W  for  each  of  Tables  2-4.  Since  the  longitude  LOP  was  also  the 
first  LOP  calculated  during  the  geolocation  technique  that  will  now  be  discussed  in 
greater  depth. 

1.  Longitude  LOP  Results 

Isolating  the  longitude  LOP  from  the  latitude  LOP  results  in  a  more  clear  picture 
of  how  accurate  the  longitude  estimation  from  the  geolocation  method  is  and  what  some 
of  the  causes  of  errors  might  be  at  a  high  level.  The  longitude  LOP  error  in  km  is  given  in 
Table  5,  with  the  negative  symbol  representing  an  estimated  position  to  the  west  of  the 
actual  position.  This  is  a  better  measure  of  accuracy  than  degrees  of  longitude,  as  the 
separation  between  them  varies  with  latitude.  The  longitude  accuracy  can  be  summarized 
as  having  a  maximum  error  of  13.697  km,  and  an  average  error  of  11.812  km  for  the 
latitude  ranges  tested.  The  errors  show  two  readily  apparent  trends,  the  first  being  that  all 
the  errors  are  from  estimated  emitter  positions  to  the  west  of  the  actual  emitter  location. 
The  next  trend  is  that  the  errors  grow  larger  in  magnitude  as  latitude  increases. 
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Table  5.  Longitude  LOP  Error  Analysis. 


Position  of  Emitter 

Longitude 
Error  (km) 

Error  In  Estimated 
Base  Frequency 
(Hz) 

Estimated  Abeam 
Time  Error  (sec) 
from  0  Hz  Doppler 
Shift  Time 

Estimated  Abeam 
Time  Error  (sec) 
from  Satellite  LEA 
Abeam  Time 

15°N 000°E 

-9.468 

16.132 

-0.18 

-1.48 

30°N  000°E 

-11.339 

32.524 

-0.45 

-3.25 

45°N  000°E 

-12.743 

53.559 

-1.06 

-5.90 

60°N  000°E 

-13.697 

79.197 

-2.54 

-10.93 

The  longitude  LOP  was  created  by  finding  the  inflection  point  on  the  received 
frequency  curve,  which  was  created  using  the  frequency  of  the  received  interference  over 
the  associated  time.  This  was  done  to  find  the  base  frequency  of  the  emitter  and  the  time 
at  which  that  base  frequency  was  received,  which  equates  to  when  the  Doppler  shift 
would  be  zero  and  the  emitter  would  be  abeam  (abeam  time).  For  a  longitude  LOP  error 
to  occur,  it  would  be  a  result  of  an  error  in  those  prior  steps. 

The  first  thing  that  was  found  in  further  analysis  was  that  there  was  an  error  in  the 
MATLAB  program’s  ability  to  calculate  the  base  frequency,  and  the  time  at  which  it 
occurred.  It  is  possible  that  an  increased  sample  rate  of  the  emitter  frequency  could  help 
this  computation,  but  that  analysis  was  not  conducted  during  this  thesis.  The  MATLAB 
calculated  base  frequency  error  is  shown  in  Table  5,  with  the  positive  denoting  that  the 
frequency  calculated  was  above  the  actual  base  frequency  of  30  GHz.  Table  5  also  shows 
the  error  in  the  time  between  where  MATLAB  calculated  the  time  of  the  inflection  point, 
as  compared  to  the  time  in  the  raw  frequency  data  that  the  Doppler  shift  became  zero. 
The  negative  for  this  column  denotes  that  the  MATLAB  calculated  inflection  point 
occurred  earlier  than  the  actual  Doppler  shift  in  the  data  became  zero.  Both  the  calculated 
frequency  and  error  in  the  abeam  time  increased  as  the  latitude  of  the  emitter  increased, 
just  as  the  position  error  does.  With  a  maximum  error  of  less  than  three  seconds  however, 
and  a  relative  rotation  velocity  between  the  satellite  and  the  emitter  of  .017  degrees/sec, 
the  maximum  longitude  error  would  be  less  than  .051  degrees.  This  meant  that  the  error 
could  not  be  fully  attributed  to  this. 
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The  emitters  were  also  all  loeated  along  the  same  meridian  (Prime  Meridian),  and 
therefore  the  abeam  time  and  the  time  at  whieh  the  Doppler  shift  beeame  zero  should  not 
have  varied  with  the  emitters’  latitudes.  In  eomparing  the  time  at  which  the  Doppler  shift 
became  zero  in  the  STK  generated  report  to  the  time  at  which  the  satellite  modeled  in 
STK  crossed  the  Prime  Meridian,  an  additional  discrepancy  was  discovered.  The  time  at 
which  the  satellite  was  at  0°E  according  to  the  latitude,  longitude,  altitude  (LLA) 
coordinates  displayed  within  STK  differed  significantly  in  each  instance  from  the  time  at 
which  the  Doppler  shift  was  zero.  In  each  case  the  time  at  which  the  Doppler  shift  was  at 
zero  was  earlier  than  the  satellite’s  longitude  was  0°E,  and  the  difference  grew 
significantly  for  the  emitters  at  higher  latitudes.  While  there  is  a  difference  in  range 
between  the  various  emitters  and  the  satellite,  the  difference  in  range  (around  4000  km)  is 
far  too  small  for  the  speed  of  light  and  transit  time  for  the  signal  (difference  of  about  .01 
seconds,  increased  from  .03  to  .04  seconds  total)  to  account  for  the  time  differences.  This 
time  difference  compounded  the  error  between  the  MATLAB  calculated  abeam  time  and 
the  zero  Doppler  shift  time,  resulting  in  the  time  differences  shown  in  the  final  column  of 
Table  5.  It  is  unclear  the  reason  for  the  differences  in  times,  but  they  are  undoubtedly  a 
significant  contributor  to  the  error  in  longitude.  At  this  time,  additional  analysis  would  be 
needed  to  understand  this  error,  and  testing  would  be  useful  to  determine  if  it  is  a  result 
only  of  the  modeling  used  in  STK  or  if  the  same  time  discrepancies  would  occur  under 
actual  conditions.  If  the  error  discussed  were  to  be  consistent  with  operational  testing 
results,  it  could  be  characterized  over  the  full  range  of  latitudes  and  then  corrected  for. 

2.  Latitude  LOP  Results 

The  error  in  km  resulting  from  the  latitude  LOP  compared  with  the  latitude  of  the 
actual  emitter  is  given  in  Table  6,  where  a  negative  sign  represents  that  the  estimated 
position  was  south  of  the  actual  position.  The  maximum  error  was  6.702  km,  with  an 
average  error  distance  of  3.776  km.  These  results  show  that  the  estimations  for  the  30°N 
emitter  were  always  the  most  accurate.  This  was  followed  by  the  45°N  emitter.  The  45°N 
and  60°N  emitter  locations  were  always  South  of  the  estimated  location  regardless  of 
where  in  the  BCA  the  emitter  was  placed.  The  15°N  and  30°N  latitude  estimates  were  not 

consistently  to  the  North  or  South  of  the  actual  emitter. 
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Table  6.  Latitude  LOP  Error  Analysis. 


Position  of 
Emitter 

Eatitude  Error  (km) 
Emitter  at  Boresight 

Eatitude  Error  (km) 
Boresight  3°S  of  Emitter 

Eatitude  Error  (km) 
Boresight  3°N  of  Emitter 

15°N  000°E 

5.680 

-6.702 

-6.702 

30°N  000°E 

-0.443 

-1.800 

0.914 

45°N  000°E 

1.569 

2.735 

2.346 

60°N  000°E 

4.954 

4.954 

6.513 

The  latitude  LOP  was  ealeulated  from  the  eomparison  of  the  reeeived  frequeney 
of  the  emitter  to  the  generated  eurves  for  simulated  emitters  at  various  assigned  latitudes. 
The  only  input  to  the  generated  curves  tied  to  the  actual  emitter  was  the  calculated  base 
frequency  based  on  the  inflection  point/abeam  time  calculation  previously  discussed.  The 
error  in  frequency,  which  was  previously  shown  in  Table  5,  is  therefore  a  potential  source 
of  error.  Another  factor  in  generating  the  curves  for  comparison  is  the  conversion  of 
geodetic  latitude  to  Cartesian  coordinates,  which  uses  an  ellipsoid  to  represent  the  Earth. 
While  this  model  is  extremely  useful,  it  is  doubtful  that  it  perfectly  represents  the  actual 
Earth,  or  the  STK  model  of  Earth  that  generated  the  frequency  data  for  the  actual  emitter 
to  be  located.  The  ellipsoid  model  could  very  well  be  the  most  reasonable  explanation  for 
the  fact  that  the  mid-latitudes  were  consistently  the  most  accurate,  if  that  is  where  the 
ellipsoid  used  most  accurately  represents  the  STK  model  (in  this  case)  or  the  real  world. 
There  exists  very  accurate  models  of  the  Earth’s  shape,  and  incorporation  of  one  into  the 
geolocation  program  could  provide  improved  accuracy  over  the  ellipsoidal  model  used. 

Eor  the  comparison  curve  generation,  in  the  current  version  of  the  program  they 
are  spaced  at  .5°  latitude  increments  and  only  created  for  3.5°S  to  3.5°N  of  the  beam’s 
boresight  latitude.  Creation  of  more  of  these  curves  and  spacing  them  closer  together 
would  provide  more  points  for  the  interpolation  step  which  estimates  the  emitter’s 
latitude.  In  the  current  method  there  are  15  points  that  a  curve  is  then  fit  to.  More  points 
could  result  in  a  more  precise  curve  fit  and  better  latitude  estimation,  although  it  does 
appear  to  be  smooth  and  already  well  characterized  by  the  curve.  The  results  did  not 
show  noticeable  trends  based  on  where  the  emitter  was  placed  inside  the  BCA,  which  is 

positive  feedback  for  the  geolocation  method  and  one  of  the  reasons  for  associating  the 
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comparison  curves  with  geodetic  latitudes.  The  eomparison  eurves  also  do  not  need  to  be 
limited  to  within  the  BCA;  they  eould  be  eomputed  as  elose  together  as  desired  from  the 
equator  to  the  northern  extent  of  the  eoverage  area.  Doing  so  may  not  inerease  aeeuraey  if 
there  are  suffieient  eurves  eomputed  within  the  BCA.  While  it  was  not  demonstrated  in 
these  results,  what  also  must  be  remembered  with  respeet  to  the  BCA  is  that  the  beam’s 
position  (north  or  south)  is  being  used  to  resolve  ambiguity  between  symmetrieal 
latitudes,  so  if  the  satellite’s  beam  was  eentered  on  the  equator  there  would  be  two 
possible  latitudes  of  the  emitter.  In  that  ease,  it  is  possible  that  a  slight  satellite  inelination 
would  result  in  a  unique  loeation,  but  that  fell  outside  the  seope  of  this  thesis. 

3.  Conclusion 

The  errors  for  the  latitude  and  longitude  LOPs  are  summarized  in  Table  7.  While 
steps  eould  be  taken  to  improve  the  aeeuraey  of  both  proeesses,  some  insight  has  been 
provided  for  where  those  errors  seem  to  originate.  Currently,  the  latitude  LOP  seems  to 
be  signifieantly  more  aceurate  than  the  longitude  LOP.  The  longitude  LOP  aeeuraey  is 
signifieantly  impaeted  by  the  timing  issue  diseussed,  whieh  may  be  a  symptom  of  the 
modeling  eondueted  and  not  persist  under  aetual  eonditions,  and  eould  be  eorreeted  for. 
Without  further  eorreetion  for  either  method,  the  results  do  prove  that  the  new  method  is 
a  viable  means  of  eondueting  geoloeation  of  EMI  emitters  using  03b  satellites. 


Table  7.  Magnitude  of  Error  Summary  for  Eatitude  and  Eongitude. 


Magnitude  of  Avg. 

Error  in  km 

Magnitude  of  Max. 

Error  in  km 

Magnitude  of  Min. 

Error  in  km 

Eatitude  EOP 

3.776 

6.702 

0.443 

Eongitude  EOP 

11.812 

13.697 

9.468 
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B.  UTILITY  OF  THE  GEOLOCATION  METHOD 

The  results  diseussed  serve  as  proof  of  concept  for  the  geolocation  technique. 
While  accuracy  could  be  improved,  in  the  current  form  the  results  would  still  prove 
useful  to  the  DOD,  or  for  that  matter,  any  user  of  03b  satellite  communications.  The 
usefulness  of  the  technique  at  the  current  accuracy  level  will  be  discussed  in  this  section. 

Two  steps  of  the  JSIR  process  are  characterization  of  the  interference  source,  and 
geolocation.75  Ugg  of  this  method  could  help  determine  whether  the  interfering  source  is 
friendly  and  interference  is  unintentional,  or  it  is  hostile  and  intentional.  In  the  case  where 
the  interfering  source  location  is  associated  with  friendly  units,  the  method  could  provide 
insight  as  to  which  unit  or  base  the  interference  is  coming  from.  In  the  event  the 
interference  is  intentional,  it  could  provide  intelligence  on  adversary  capabilities  or 
support  networks  used  to  acquire  the  interference  equipment.  If  in  a  complicated  area  of 
operations  with  multiple  actors,  the  geolocation  method  could  place  the  location  inside 
territory  controlled  by  one  of  them,  again  providing  useful  knowledge.  In  the  event  the 
interference  appears  to  be  the  result  of  a  state  actor,  it  could  enable  diplomatic  actions  to 
be  taken  to  cease  the  interference. 

If  the  interference  cannot  be  stopped  at  the  source,  a  simple  way  interference  can 
be  resolved  is  to  re-point  the  BCA  containing  the  interference  such  that  the  interfering 
source  is  inside  a  null  of  the  antenna.  The  narrow  beam  width  is  one  of  the  advantages  of 
Ka-band,  which  can  be  used  for  exactly  this  purpose.  While  there  is  significant  gain  out 
to  the  -6  dB  point  in  the  03b  payload  antenna  receive  pattern,  beyond  that  the  gain  drops 
off  sharply,  as  was  shown  in  Figure  1  and  Figure  4.  Moving  the  beam  slightly  can  place 
the  interfering  source  into  this  disadvantageous  null  in  the  pattern  while  maintaining  the 
user  in  the  higher  gain  of  the  beam  pattern. 

While  the  level  of  accuracy  achieved  by  this  method  is  useful  on  its  own,  it  can 
also  be  enhanced  by  using  other  knowledge  and  situations  that  may  be  applicable.  An 
example  is  that  if  the  location  area  includes  water,  it  may  be  more  likely  that  the 
interfering  source  is  located  on  the  land  nearby.  Another  situation  would  be  that  if  there 

“Joint  Spectrum  Interference  Resolution  (JSIR)  Procedures,”  A-6  (see  chap.  1,  n.  14). 


58 


have  been  interfering  sources  that  were  previously  located  by  more  accurate  means,  a 
location  area  including  that  source  location  could  possibly  be  assumed  to  be  a  repeat 
offense  by  the  same  source.  In  a  remote  area,  there  may  be  little  development  which 
could  hide  an  interfering  antenna/source.  Finally,  the  United  States  or  its  allies  may  have 
assets  capable  of  providing  a  more  accurate  location  of  the  interfering  source.  The  initial 
estimate  provided  by  the  geolocation  method  developed  may  allow  these  other  assets  to 
focus  their  search  area  so  that  they  may  be  used  most  efficiently,  ensuring  they  are 
available  for  other  tasking  as  soon  as  possible.  If  the  worst  case  geolocation  is 
conservatively  considered  to  have  a  radius  for  a  circular  area  of  probability  (AOP)  of  20 
km,  and  the  BCA  where  the  emitter  can  be  assumed  to  be  located  has  a  radius  of  350  km, 
the  area  needed  to  be  searched  by  more  capable  assets  has  been  reduced  by  99.7% 
through  use  of  this  method,  greatly  increasing  their  efficiency  for  what  is  likely  a  low 
density  asset.  While  any  of  these  supplementary  procedures  are  pending  or  taking  place, 
03b  can  easily  restore  communications  by  beam  steering,  as  discussed  in  the  previous 
paragraph. 

C.  FUTURE  WORK 

While  this  thesis  did  achieve  a  proof  of  concept  for  a  new  geolocation  method  in 
Ka-band  using  03b  satellites,  there  is  still  opportunity  for  further  refinement.  The  first  of 
these  opportunities  is  to  conduct  a  sensitivity  analysis  that  would  identify  the  largest 
causes  of  errors.  Among  known  errors  and  areas  for  refinement  is  the  ability  of  the 
program  to  calculate  the  base  frequency  of  the  transmitter  from  the  second  derivative  of 
the  received  frequency  curve.  This  error,  with  the  abeam  time  calculation,  results  in  a 
portion  of  the  error  in  the  longitude  LOP  calculation.  A  much  larger  portion  of  the  error 
was  traced  back  to  STK  and  a  time  discrepancy  between  when  the  satellite  was  abeam  the 
emitter  and  the  time  at  which  the  Doppler  shift  went  to  zero,  with  no  apparent  cause, 
warranting  further  research.  For  the  latitude  LOP,  creating  additional  comparison  curves 
and  optimizing  the  curve  fit  for  the  interpolation  curve  resulting  from  the  added  points 
may  result  in  improved  accuracy.  Along  with  error  resolution/accuracy  enhancement 
should  be  adding  the  last  increment  of  realism  to  the  modeling,  which  would  involve 

using  the  03b  satellite  ephemeris  data  instead  of  the  simplified  orbit  at  0°  inclination  and 
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no  eccentricity.  In  addition  to  this,  this  thesis  has  assumed  that  the  carrier  frequency  of 
the  interfering  signal  can  be  monitored  to  the  same  level  of  precision  simulated  by  the 
STK  scenario.  While  03b  Networks  does  use  carrier  monitoring  software  at  gateway 
locations,  further  research  would  be  needed  to  see  if  this  currently  provides  the  data 
necessary  to  implement  use  of  this  geolocation  method.^^  If  the  current  system  does  not 
provide  the  data  necessary,  then  signal  processing  may  be  necessary  to  obtain  the  data 
required  for  practical  use.  Actual  testing  of  the  method  would  also  be  a  necessary  step 
before  implementation  could  be  considered.  Testing  could  then  enable  a  thorough  error 
analysis  under  realistic  conditions  and  not  under  the  constraints  of  modeling,  and 
appropriate  calibration  applied  to  optimize  use  of  the  method. 

While  not  directly  related  to  this  method,  a  different  approach  to  using  03b 
satellites  to  perform  the  geolocation  was  briefly  discussed  with  thesis  advisors.  This 
would  involve  treating  the  geolocation  as  a  general  estimation  problem,  which  could  then 
be  solved  using  an  iterative  method  (such  as  the  Newton-Raphson  method).  Due  to  time 
constraints  this  approach  was  not  pursued,  but  it  could  be  an  alternative  to  the  developed 
method  and  provide  the  same  benefits  of  availability  and  low  cost  with  a  useful  level  of 
accuracy. 

D.  CONCLUSION 

For  the  DOD  this  geolocation  method  has  great  potential.  Likely  the  two  greatest 
benefits  to  the  DOD  are  cost  and  availability.  This  method  can  be  executed  using 
MATLAB,  Microsoft  Excel,  and  an  inexpensive  laptop,  all  of  which  could  be  made 
readily  available  by  units  in  the  DOD.  These  can  be  reasonably  assumed  to  be  lower  cost 
than  a  dedicated  signals  intelligence  (SIGINT)  or  electronic  warfare  (EW)  asset 
providing  geolocation  capability.  The  cost  of  a  laptop  with  the  identified  programs  might 
run  hundreds  of  dollars,  whereas  a  single  plane  could  cost  millions  of  dollars,  and  for  a 
satellite,  hundreds  of  millions  of  dollars.  Eurthermore,  there  is  no  risk  involved,  as  no 
asset  needs  to  be  deployed  to  conduct  the  geolocation  unless  additional  accuracy  is 
desired.  The  only  input  necessary  for  the  method  to  work  is  a  received  frequency  and 

Ken  Mentasti,  email  message  to  author,  February  1,  2016. 
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associated  time,  which  would  already  be  monitored  by  03b  Networks  or/and  the 
eustomer.  Use  of  03b  for  eommunieations  would  presumably  already  be  paid  for,  as  the 
DOD  would  not  be  concerned  if  an  03b  satellite  experieneed  jamming  unless  the  DOD 
was  the  customer.  03b  Networks  may  even  desire  to  implement  the  geoloeation  method 
in  order  to  provide  customers  better  service  by  mitigation  of  interference  using  beam 
steering.  If  the  method  was  implemented  by  the  eompany  instead  of  the  DOD,  then  the 
ability  to  obtain  the  geolocation  from  them  to  gain  any  potential  intelligence  benefit  from 
it  eould  surely  be  worked  out.  The  accuraey  in  the  method’s  current  form  is  enough  to  be 
useful,  but  may  not  be  an  improvement  over  existing  abilities.  While  that  is  true,  an  asset 
that  is  capable  of  greater  aeeuracy  but  not  available  to  the  user  at  the  time  interference 
oecurs,  is  not  useful  at  all.  The  new  geolocation  method  does  require  for  the  satellite  to 
pass  abeam  the  emitter  before  the  geoloeation  ean  be  obtained,  which  does  take  time  but 
is  likely  still  faster  than  a  more  capable  (and  typieally  expensive  and  scarce)  asset  eould 
be  tasked.  The  timeliness  of  this  geoloeation  would  also  not  depend  on  the  unit’s  priority 
for  the  tasking  of  the  more  capable  asset. 

This  geoloeation  method  eould  beeome  an  organie,  low  cost,  100%  available 
capability  to  any  unit  with  the  ability  to  run  the  reeeived  frequency  and  time  data  of  the 
interference  souree  through  the  MATLAB  program.  The  timeliness  of  obtaining  an 
estimate  for  the  emitter  location  would  only  depend  on  how  long  it  took  for  the  satellite 
to  transit  abeam  the  emitter.  The  aceuraey  is  currently  enough  to  mitigate  the  interferenee 
and  possibly  provide  useful  intelligence,  with  the  potential  for  the  accuracy  to  be 
inereased  beyond  the  current  state.  These  are  substantial  benefits  to  the  DOD  whieh  are 
eurrently  not  provided  by  other  means  and  are  worth  pursuing — partieularly  in  light  of 
the  erueial  role  SATCOM  plays  for  the  DOD  and  the  growing  use  of  Ka-band 
frequencies  in  that  role. 


61 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


62 


APPENDIX  A.  SAMPLE  DATA  EOR  USE  IN  MATLAB  PROGRAM 


Table  8.  STK  Generated  Data  for  Emitter  at  30°N  000°E. 


Time  (UTCG) 

Reev  Interferenee  Elapsed  Time 

Revd.  Erequeney  (Hz) 

29  Mar  2016  20:17:20.000 

347 

30000010947.08560 

29  Mar  2016  20:17:21.000 

348 

30000010875.94140 

29  Mar  2016  20:17:22.000 

349 

30000010804.79390 

29  Mar  2016  20:17:23.000 

350 

30000010733.64300 

29  Mar  2016  20:17:24.000 

351 

30000010662.48870 

29  Mar  2016  20:17:25.000 

352 

30000010591.33120 

29  Mar  2016  20:17:26.000 

353 

30000010520.17040 

29  Mar  2016  20:17:27.000 

354 

30000010449.00630 

29  Mar  2016  20:17:28.000 

355 

30000010377.83910 

29  Mar  2016  20:17:29.000 

356 

30000010306.66860 

29  Mar  2016  20:17:30.000 

357 

30000010235.49490 

29  Mar  2016  20:17:31.000 

358 

30000010164.31810 

29  Mar  2016  20:17:32.000 

359 

30000010093.13810 

29  Mar  2016  20:17:33.000 

360 

30000010021.95500 

29  Mar  2016  20:17:34.000 

361 

30000009950.76880 

29  Mar  2016  20:17:35.000 

362 

30000009879.57950 

29  Mar  2016  20:17:36.000 

363 

30000009808.38720 

29  Mar  2016  20:17:37.000 

364 

30000009737.19180 

29  Mar  2016  20:17:38.000 

365 

30000009665.99340 

29  Mar  2016  20:17:39.000 

366 

30000009594.79210 

29  Mar  2016  20:17:40.000 

367 

30000009523.58770 

29  Mar  2016  20:17:41.000 

368 

30000009452.38050 

29  Mar  2016  20:17:42.000 

369 

30000009381.17030 

29  Mar  2016  20:17:43.000 

370 

30000009309.95720 

29  Mar  2016  20:17:44.000 

371 

30000009238.74120 

29  Mar  2016  20:17:45.000 

372 

30000009167.52240 

29  Mar  2016  20:17:46.000 

373 

30000009096.30070 

29  Mar  2016  20:17:47.000 

374 

30000009025.07620 

29  Mar  2016  20:17:48.000 

375 

30000008953.84890 

29  Mar  2016  20:17:49.000 

376 

30000008882.61890 

29  Mar  2016  20:17:50.000 

377 

30000008811.38610 

29  Mar  2016  20:17:51.000 

378 

30000008740.15060 

29  Mar  2016  20:17:52.000 

379 

30000008668.91240 

29  Mar  2016  20:17:53.000 

380 

30000008597.67150 
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29  Mar  2016  20:17:54.000 

381 

30000008526.42790 

29  Mar  2016  20:17:55.000 

382 

30000008455.18170 

29  Mar  2016  20:17:56.000 

383 

30000008383.93290 

29  Mar  2016  20:17:57.000 

384 

30000008312.68150 

29  Mar  2016  20:17:58.000 

385 

30000008241.42760 

29  Mar  2016  20:17:59.000 

386 

30000008170.17110 

29  Mar  2016  20:18:00.000 

387 

30000008098.91200 

29  Mar  2016  20:18:01.000 

388 

30000008027.65050 

29  Mar  2016  20:18:02.000 

389 

30000007956.38640 

29  Mar  2016  20:18:03.000 

390 

30000007885.11990 

29  Mar  2016  20:18:04.000 

391 

30000007813.85100 

29  Mar  2016  20:18:05.000 

392 

30000007742.57970 

29  Mar  2016  20:18:06.000 

393 

30000007671.30590 

29  Mar  2016  20:18:07.000 

394 

30000007600.02980 

29  Mar  2016  20:18:08.000 

395 

30000007528.75140 

29  Mar  2016  20:18:09.000 

396 

30000007457.47060 

29  Mar  2016  20:18:10.000 

397 

30000007386.18750 

29  Mar  2016  20:18:11.000 

398 

30000007314.90210 

29  Mar  2016  20:18:12.000 

399 

30000007243.61450 

29  Mar  2016  20:18:13.000 

400 

30000007172.32460 

29  Mar  2016  20:18:14.000 

401 

30000007101.03250 

29  Mar  2016  20:18:15.000 

402 

30000007029.73820 

29  Mar  2016  20:18:16.000 

403 

30000006958.44170 

29  Mar  2016  20:18:17.000 

404 

30000006887.14310 

29  Mar  2016  20:18:18.000 

405 

30000006815.84230 

29  Mar  2016  20:18:19.000 

406 

30000006744.53950 

29  Mar  2016  20:18:20.000 

407 

30000006673.23460 

29  Mar  2016  20:18:21.000 

408 

30000006601.92760 

29  Mar  2016  20:18:22.000 

409 

30000006530.61850 

29  Mar  2016  20:18:23.000 

410 

30000006459.30740 

29  Mar  2016  20:18:24.000 

411 

30000006387.99440 

29  Mar  2016  20:18:25.000 

412 

30000006316.67930 

29  Mar  2016  20:18:26.000 

413 

30000006245.36240 

29  Mar  2016  20:18:27.000 

414 

30000006174.04350 

29  Mar  2016  20:18:28.000 

415 

30000006102.72260 

29  Mar  2016  20:18:29.000 

416 

30000006031.39990 

29  Mar  2016  20:18:30.000 

417 

30000005960.07540 

29  Mar  2016  20:18:31.000 

418 

30000005888.74890 

29  Mar  2016  20:18:32.000 

419 

30000005817.42070 

29  Mar  2016  20:18:33.000 

420 

30000005746.09070 
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29  Mar  2016  20:18:34.000 

421 

30000005674.75880 

29  Mar  2016  20:18:35.000 

422 

30000005603.42530 

29  Mar  2016  20:18:36.000 

423 

30000005532.09000 

29  Mar  2016  20:18:37.000 

424 

30000005460.75300 

29  Mar  2016  20:18:38.000 

425 

30000005389.41430 

29  Mar  2016  20:18:39.000 

426 

30000005318.07390 

29  Mar  2016  20:18:40.000 

427 

30000005246.73190 

29  Mar  2016  20:18:41.000 

428 

30000005175.38830 

29  Mar  2016  20:18:42.000 

429 

30000005104.04310 

29  Mar  2016  20:18:43.000 

430 

30000005032.69630 

29  Mar  2016  20:18:44.000 

431 

30000004961.34800 

29  Mar  2016  20:18:45.000 

432 

30000004889.99810 

29  Mar  2016  20:18:46.000 

433 

30000004818.64670 

29  Mar  2016  20:18:47.000 

434 

30000004747.29380 

29  Mar  2016  20:18:48.000 

435 

30000004675.93950 

29  Mar  2016  20:18:49.000 

436 

30000004604.58370 

29  Mar  2016  20:18:50.000 

437 

30000004533.22650 

29  Mar  2016  20:18:51.000 

438 

30000004461.86790 

29  Mar  2016  20:18:52.000 

439 

30000004390.50790 

29  Mar  2016  20:18:53.000 

440 

30000004319.14650 

29  Mar  2016  20:18:54.000 

441 

30000004247.78390 

29  Mar  2016  20:18:55.000 

442 

30000004176.41990 

29  Mar  2016  20:18:56.000 

443 

30000004105.05460 

29  Mar  2016  20:18:57.000 

444 

30000004033.68810 

29  Mar  2016  20:18:58.000 

445 

30000003962.32030 

29  Mar  2016  20:18:59.000 

30000003890.95130 

29  Mar  2016  20:19:00.000 

447 

30000003819.58110 

29  Mar  2016  20:19:01.000 

448 

30000003748.20970 

29  Mar  2016  20:19:02.000 

449 

30000003676.83720 

29  Mar  2016  20:19:03.000 

450 

30000003605.46350 

29  Mar  2016  20:19:04.000 

451 

30000003534.08880 

29  Mar  2016  20:19:05.000 

452 

30000003462.71290 

29  Mar  2016  20:19:06.000 

453 

30000003391.33600 

29  Mar  2016  20:19:07.000 

454 

30000003319.95800 

29  Mar  2016  20:19:08.000 

455 

30000003248.57900 

29  Mar  2016  20:19:09.000 

456 

30000003177.19900 

29  Mar  2016  20:19:10.000 

457 

30000003105.81800 

29  Mar  2016  20:19:11.000 

458 

30000003034.43610 

29  Mar  2016  20:19:12.000 

459 

30000002963.05320 

29  Mar  2016  20:19:13.000 

460 

30000002891.66940 
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29  Mar  2016  20:19:14.000 

461 

30000002820.28480 

29  Mar  2016  20:19:15.000 

462 

30000002748.89920 

29  Mar  2016  20:19:16.000 

463 

30000002677.51280 

29  Mar  2016  20:19:17.000 

464 

30000002606.12560 

29  Mar  2016  20:19:18.000 

465 

30000002534.73760 

29  Mar  2016  20:19:19.000 

466 

30000002463.34890 

29  Mar  2016  20:19:20.000 

467 

30000002391.95930 

29  Mar  2016  20:19:21.000 

468 

30000002320.56910 

29  Mar  2016  20:19:22.000 

469 

30000002249.17810 

29  Mar  2016  20:19:23.000 

470 

30000002177.78640 

29  Mar  2016  20:19:24.000 

471 

30000002106.39410 

29  Mar  2016  20:19:25.000 

472 

30000002035.00110 

29  Mar  2016  20:19:26.000 

473 

30000001963.60750 

29  Mar  2016  20:19:27.000 

474 

30000001892.21330 

29  Mar  2016  20:19:28.000 

475 

30000001820.81850 

29  Mar  2016  20:19:29.000 

476 

30000001749.42310 

29  Mar  2016  20:19:30.000 

477 

30000001678.02730 

29  Mar  2016  20:19:31.000 

478 

30000001606.63090 

29  Mar  2016  20:19:32.000 

479 

30000001535.23400 

29  Mar  2016  20:19:33.000 

480 

30000001463.83670 

29  Mar  2016  20:19:34.000 

481 

30000001392.43890 

29  Mar  2016  20:19:35.000 

482 

30000001321.04070 

29  Mar  2016  20:19:36.000 

483 

30000001249.64200 

29  Mar  2016  20:19:37.000 

484 

30000001178.24310 

29  Mar  2016  20:19:38.000 

485 

30000001106.84370 

29  Mar  2016  20:19:39.000 

486 

30000001035.44400 

29  Mar  2016  20:19:40.000 

487 

30000000964.04400 

29  Mar  2016  20:19:41.000 

488 

30000000892.64380 

29  Mar  2016  20:19:42.000 

489 

30000000821.24320 

29  Mar  2016  20:19:43.000 

490 

30000000749.84240 

29  Mar  2016  20:19:44.000 

491 

30000000678.44140 

29  Mar  2016  20:19:45.000 

492 

30000000607.04020 

29  Mar  2016  20:19:46.000 

493 

30000000535.63880 

29  Mar  2016  20:19:47.000 

494 

30000000464.23720 

29  Mar  2016  20:19:48.000 

495 

30000000392.83550 

29  Mar  2016  20:19:49.000 

496 

30000000321.43370 

29  Mar  2016  20:19:50.000 

497 

30000000250.03180 

29  Mar  2016  20:19:51.000 

498 

30000000178.62990 

29  Mar  2016  20:19:52.000 

499 

30000000107.22790 

29  Mar  2016  20:19:53.000 

500 

30000000035.82580 
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29  Mar  2016  20:19:54.000 

501 

29999999964.42380 

29  Mar  2016  20:19:55.000 

502 

29999999893.02180 

29  Mar  2016  20:19:56.000 

503 

29999999821.61980 

29  Mar  2016  20:19:57.000 

504 

29999999750.21790 

29  Mar  2016  20:19:58.000 

505 

29999999678.81610 

29  Mar  2016  20:19:59.000 

506 

29999999607.41440 

29  Mar  2016  20:20:00.000 

507 

29999999536.01290 

29  Mar  2016  20:20:01.000 

508 

29999999464.61140 

29  Mar  2016  20:20:02.000 

509 

29999999393.21020 

29  Mar  2016  20:20:03.000 

510 

29999999321.80920 

29  Mar  2016  20:20:04.000 

511 

29999999250.40830 

29  Mar  2016  20:20:05.000 

512 

29999999179.00780 

29  Mar  2016  20:20:06.000 

513 

29999999107.60750 

29  Mar  2016  20:20:07.000 

514 

29999999036.20740 

29  Mar  2016  20:20:08.000 

515 

29999998964.80770 

29  Mar  2016  20:20:09.000 

516 

29999998893.40830 

29  Mar  2016  20:20:10.000 

517 

29999998822.00930 

29  Mar  2016  20:20:1 1.000 

518 

29999998750.61070 

29  Mar  2016  20:20:12.000 

519 

29999998679.21240 

29  Mar  2016  20:20:13.000 

520 

29999998607.81460 

29  Mar  2016  20:20:14.000 

521 

29999998536.41720 

29  Mar  2016  20:20:15.000 

522 

29999998465.02030 

29  Mar  2016  20:20:16.000 

523 

29999998393.62380 

29  Mar  2016  20:20:17.000 

524 

29999998322.22790 

29  Mar  2016  20:20:18.000 

525 

29999998250.83250 

29  Mar  2016  20:20:19.000 

526 

29999998179.43770 

29  Mar  2016  20:20:20.000 

527 

29999998108.04340 

29  Mar  2016  20:20:21.000 

528 

29999998036.64980 

29  Mar  2016  20:20:22.000 

529 

29999997965.25670 

29  Mar  2016  20:20:23.000 

530 

29999997893.86430 

29  Mar  2016  20:20:24.000 

531 

29999997822.47260 

29  Mar  2016  20:20:25.000 

532 

29999997751.08150 

29  Mar  2016  20:20:26.000 

533 

29999997679.69120 

29  Mar  2016  20:20:27.000 

534 

29999997608.30160 

29  Mar  2016  20:20:28.000 

535 

29999997536.91270 

29  Mar  2016  20:20:29.000 

536 

29999997465.52470 

29  Mar  2016  20:20:30.000 

537 

29999997394.13740 

29  Mar  2016  20:20:31.000 

538 

29999997322.75090 

29  Mar  2016  20:20:32.000 

539 

29999997251.36530 

29  Mar  2016  20:20:33.000 

540 

29999997179.98060 
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29  Mar  2016  20:20:34.000 

541 

29999997108.59670 

29  Mar  2016  20:20:35.000 

542 

29999997037.21370 

29  Mar  2016  20:20:36.000 

543 

29999996965.83170 

29  Mar  2016  20:20:37.000 

544 

29999996894.45060 

29  Mar  2016  20:20:38.000 

545 

29999996823.07050 

29  Mar  2016  20:20:39.000 

546 

29999996751.69140 

29  Mar  2016  20:20:40.000 

547 

29999996680.31340 

29  Mar  2016  20:20:41.000 

548 

29999996608.93630 

29  Mar  2016  20:20:42.000 

549 

29999996537.56040 

29  Mar  2016  20:20:43.000 

550 

29999996466.18550 

29  Mar  2016  20:20:44.000 

551 

29999996394.81170 

29  Mar  2016  20:20:45.000 

552 

29999996323.43910 

29  Mar  2016  20:20:46.000 

553 

29999996252.06760 

29  Mar  2016  20:20:47.000 

554 

29999996180.69730 

29  Mar  2016  20:20:48.000 

555 

29999996109.32810 

29  Mar  2016  20:20:49.000 

556 

29999996037.96030 

29  Mar  2016  20:20:50.000 

557 

29999995966.59360 

29  Mar  2016  20:20:51.000 

558 

29999995895.22820 

29  Mar  2016  20:20:52.000 

559 

29999995823.86410 

29  Mar  2016  20:20:53.000 

560 

29999995752.50130 

29  Mar  2016  20:20:54.000 

561 

29999995681.13980 

29  Mar  2016  20:20:55.000 

562 

29999995609.77970 

29  Mar  2016  20:20:56.000 

563 

29999995538.42100 

29  Mar  2016  20:20:57.000 

564 

29999995467.06370 

29  Mar  2016  20:20:58.000 

565 

29999995395.70770 

29  Mar  2016  20:20:59.000 

566 

29999995324.35330 

29  Mar  2016  20:21:00.000 

567 

29999995253.00020 

29  Mar  2016  20:21:01.000 

568 

29999995181.64870 

29  Mar  2016  20:21:02.000 

569 

29999995110.29870 

29  Mar  2016  20:21:03.000 

570 

29999995038.95020 

29  Mar  2016  20:21:04.000 

571 

29999994967.60320 

29  Mar  2016  20:21:05.000 

572 

29999994896.25790 

29  Mar  2016  20:21:06.000 

573 

29999994824.91410 

29  Mar  2016  20:21:07.000 

574 

29999994753.57190 

29  Mar  2016  20:21:08.000 

575 

29999994682.23140 

29  Mar  2016  20:21:09.000 

576 

29999994610.89260 

29  Mar  2016  20:21:10.000 

577 

29999994539.55540 

29  Mar  2016  20:21:11.000 

578 

29999994468.22000 

29  Mar  2016  20:21:12.000 

579 

29999994396.88620 

29  Mar  2016  20:21:13.000 

580 

29999994325.55420 
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29  Mar  2016  20:21:14.000 

581 

29999994254.22400 

29  Mar  2016  20:21:15.000 

582 

29999994182.89560 

29  Mar  2016  20:21:16.000 

583 

29999994111.56910 

29  Mar  2016  20:21:17.000 

584 

29999994040.24430 

29  Mar  2016  20:21:18.000 

585 

29999993968.92140 

29  Mar  2016  20:21:19.000 

586 

29999993897.60040 

29  Mar  2016  20:21:20.000 

587 

29999993826.28130 

29  Mar  2016  20:21:21.000 

588 

29999993754.96410 

29  Mar  2016  20:21:22.000 

589 

29999993683.64890 

29  Mar  2016  20:21:23.000 

590 

29999993612.33570 

29  Mar  2016  20:21:24.000 

591 

29999993541.02440 

29  Mar  2016  20:21:25.000 

592 

29999993469.71520 

29  Mar  2016  20:21:26.000 

593 

29999993398.40800 

29  Mar  2016  20:21:27.000 

594 

29999993327.10280 

29  Mar  2016  20:21:28.000 

595 

29999993255.79980 

29  Mar  2016  20:21:29.000 

596 

29999993184.49880 

29  Mar  2016  20:21:30.000 

597 

29999993113.20000 

29  Mar  2016  20:21:31.000 

598 

29999993041.90340 

29  Mar  2016  20:21:32.000 

599 

29999992970.60880 

29  Mar  2016  20:21:33.000 

600 

29999992899.31650 

29  Mar  2016  20:21:34.000 

601 

29999992828.02640 

29  Mar  2016  20:21:35.000 

602 

29999992756.73860 

29  Mar  2016  20:21:36.000 

603 

29999992685.45300 

29  Mar  2016  20:21:37.000 

604 

29999992614.16970 

29  Mar  2016  20:21:38.000 

605 

29999992542.88870 

29  Mar  2016  20:21:39.000 

606 

29999992471.61000 

29  Mar  2016  20:21:40.000 

607 

29999992400.33360 

29  Mar  2016  20:21:41.000 

608 

29999992329.05970 

29  Mar  2016  20:21:42.000 

609 

29999992257.78810 

29  Mar  2016  20:21:43.000 

610 

29999992186.51900 

29  Mar  2016  20:21:44.000 

611 

29999992115.25220 

29  Mar  2016  20:21:45.000 

612 

29999992043.98800 

29  Mar  2016  20:21:46.000 

613 

29999991972.72620 

29  Mar  2016  20:21:47.000 

614 

29999991901.46690 

29  Mar  2016  20:21:48.000 

615 

29999991830.21010 

29  Mar  2016  20:21:49.000 

616 

29999991758.95590 

29  Mar  2016  20:21:50.000 

617 

29999991687.70430 

29  Mar  2016  20:21:51.000 

618 

29999991616.45520 

29  Mar  2016  20:21:52.000 

619 

29999991545.20880 

29  Mar  2016  20:21:53.000 

620 

29999991473.96500 
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29  Mar  2016  20:21:54.000 

621 

29999991402.72380 

29  Mar  2016  20:21:55.000 

622 

29999991331.48540 

29  Mar  2016  20:21:56.000 

623 

29999991260.24960 

29  Mar  2016  20:21:57.000 

624 

29999991189.01660 

29  Mar  2016  20:21:58.000 

625 

29999991117.78620 

29  Mar  2016  20:21:59.000 

626 

29999991046.55870 

29  Mar  2016  20:22:00.000 

627 

29999990975.33400 

29  Mar  2016  20:22:01.000 

628 

29999990904.11200 

29  Mar  2016  20:22:02.000 

629 

29999990832.89290 

29  Mar  2016  20:22:03.000 

630 

29999990761.67670 

29  Mar  2016  20:22:04.000 

631 

29999990690.46330 

29  Mar  2016  20:22:05.000 

632 

29999990619.25280 

29  Mar  2016  20:22:06.000 

633 

29999990548.04530 

29  Mar  2016  20:22:07.000 

634 

29999990476.84070 

29  Mar  2016  20:22:08.000 

635 

29999990405.63900 

29  Mar  2016  20:22:09.000 

636 

29999990334.44040 

29  Mar  2016  20:22:10.000 

637 

29999990263.24470 

29  Mar  2016  20:22:1 1.000 

638 

29999990192.05210 

29  Mar  2016  20:22:12.000 

639 

29999990120.86250 

29  Mar  2016  20:22:13.000 

640 

29999990049.67600 

29  Mar  2016  20:22:14.000 

641 

29999989978.49260 

29  Mar  2016  20:22:15.000 

642 

29999989907.31240 

29  Mar  2016  20:22:16.000 

643 

29999989836.13520 

29  Mar  2016  20:22:17.000 

644 

29999989764.96120 

29  Mar  2016  20:22:18.000 

645 

29999989693.79040 

29  Mar  2016  20:22:19.000 

646 

29999989622.62290 

29  Mar  2016  20:22:20.000 

647 

29999989551.45850 
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APPENDIX  B.  MATLAB  GEOLOCATION  PROGRAM  CODE 


%Program  for  Geolocation  of  EMI  Emitters  by  03b  Satellites  -  Connolly 
%Sec  1 .  Inputs 

filename  =  'ThesisDopShift.xlsx';  %Eile  name  for  Excel  stored  STK  generated  data 
sheet  =  'ElatSON';  %Sheet  Name  for  Excel  Data 
ETime=xlsread(filename, sheet, 'B3;B303');  %Read  data  from  Excel 
YJam=(xlsread(fdename, sheet, 'C3:C303'));  %  Read  data  from  Excel 
SatInitTime=4293;  %STK  Scenario  Time  at  Start  of  EMI  recv. 

%Beam  Center  (deg) 
lat=30; 

%Eor  Error  Calculator 
ActEat=30; 

ActEong=0; 

%Sat  Eongitude  at  Scenario  Start 
SatInitEong=-80; 

%Time  Step  for  Baseline  Rdot  Calcs 
dt=.01; 

%Sec  2.  Constants  and  Properties 
%mu 

mu  =  3.986004el4; 

%Earth's  radii  (m) 

Re=6378136; 

%Oblate  Earth 
e=081819301; 
ellipsoid=[Re  e]; 

%Sidereal  day  (s) 
sd=86 164.1; 

%Speed  of  light  (m/s) 
c  =  3e8; 

%Satellite  Altitude  (m) 

Z  =  8063e3; 

%Semimaior  Axis  (m) 
a=Re+Z; 

%Sat.  Orbital  Period  (s) 

T  =  2*pi*(sqrt(a.^3/mu)); 

%Relative  Angular  Velocity  of  Sat  to  Earth  (rad/s) 

RotV=(2*pi)/(T)-(2*pi)/(sd); 
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%Sec  3.  Jammer  abeam  time  ealeulation,  eentering,  plotting 

RawDat=polyfit(ETime,  YJam,  3); 

Dfreq=polyder(RawDat); 

DDfreq=polyder(Dfreq); 

AbeamTime=roots(DDfreq); 

FnotJam=polyval(RawDat,  AbeamTime); 

XRaw=-ETime; 

TJam=XRaw+AbeamTime;  %Center  data  at  t=0 
maxT  ime=max(T  J  am) ; 
mini  ime=min(T  Jam); 

ReeEreqJam=polyfit(TJam,YJam,3); 

figure 

plot  (TJam,  YJam,'— k'); 

title('Reoeived  Ereq  (Hz)  Relative  to  Satellite  Eoeation') 
xlabel('Time  (s)  Relative  to  Abeam  Point') 
ylabel('Reoeived  Erequeney  (Hz)') 
hold  on 

Y  JamFit=poly  val(ReoE  req  J  am,  T  J  am) ; 

Y  J  amResid= Y  J  am- Y  J  amF  it; 

SSresid=sum(YJamResidA2); 

SStotal=(length(YJam)-l)*var(YJam); 

rsq=l-SSresid/SStotal; 

adj_rsq=l-SSresid/SStotal*(length(YJam)-l)/(length(YJam)-... 
length(ReoF  req  J  am)); 

XTraekJam=polyval(ReeFreqJam,maxTime); 

=lc  He  *  =[c  =lc  He  *  He  =lc  =[C  *  =1=  *=[=*=1=  ***=[=  ***=1=  ***=[==1=  ****=[=*=[=  ***=1=  *****=[=*  **=[=*  ****=1=  ***=[==1=  ***=1=  ***=1=  * 

%See.  4.  Eongitude  LOP  ealeulation 

LongEst=SatInitLong+(RotV*(SatInitTime+AbeamTime)*  1 80/pi) 

%Seo  5.  Ree.  Ereq  Curves  of  Sim  Emitters  for  Comparison/Latitude  LOP 

Tbase=(minTime:dt:maxTime); 

T  dbase=(minT  ime;  dt:  (maxT  ime+dt)); 

Angle=RotV*Tdbase; 

%Seo.  5.1.  Baseline  Curve  3.5  deg  Below  Beam  Center 
latl=(lat-3.5)*pi/180;  %deg  to  rad 
Nphi  1  =Re/sqrt(  1  -e.^2  *  (sin(lat  1 ))  .^2); 
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JamPosXl=a-Nphil*cos(latl)*cos(  Angle); 

JamPosYl=Nphil  *cos(latl)*sin(Angle); 

JamPosZ  1 =(  1 -e  A2)  *Nphi  1  *  sin(lat  1 ); 
JamPosl=sqrt(JamPosXl  A2+JamPosY  1  A2+JamPosZl  A2); 
Rdot  1  =(diff(  JamPos  1 ))/ dt; 

Yl=(Rdotl/c)*FnotJam+FnotJam; 
plot  (Tbase,  Yl); 
hold  on 

RecFreql=polyfit(Tbase,Yl,3); 

XTrack  1  =poly  val(RecFreq  1  ,maxT  ime); 
YJaniFitl=polyval(RecFreql, Tbase); 

Y  JamResid  1  =Y  1  -  Y  JamFit  1 ; 

SSresidl=sum(YJamResidl  A2); 

SStotal  1  =(length(Y  1 )- 1  )*  var(  Y 1 ); 
rsql=l-SSresidl./SStotall; 

adj_rsq  1 = 1  -S  Sresid  1  ./S  Stotal  1  *(length(  Y 1 )- 1 )  ./(length(Y  1 )- . . . 
length(ReeF  req  1 )) ; 

%See.  5.2.  Baseline  Curve  3  deg  Below  Beam  Center 
lat2=(lat-3)*pi/180;  %deg  to  rad 
Nphi2=Re/sqrt(  1  -e.^2  *  (sin(lat2))  .^2); 
JamPosX2=a-Nphi2*oos(lat2)*oos(  Angle); 
JamPosY2=Nphi2*eos(lat2)*sin(  Angle); 
JamPosZ2=(l-e.^2)*Nphi2*sin(lat2); 
JamPos2=sqrt(JamPosX2.^2+JamPosY2.^2+JamPosZ2.^2); 
Rdot2=(diff(JamPos2))/dt; 

Y  2=(Rdot2/ c)  *  Fnot  J  am+F  not  J  am; 
plot  (Tbase,  Y2); 

hold  on 

RecFreq2=polyfit(Tbase,Y2,3); 
XTrack2=polyval(ReoFreq2,maxTime); 
YJamFit2=polyval(ReoFreq2, Tbase); 

Y  JamResid2=Y2- Y  JamFit2; 

S  Sresid2=sum(  Y  J  amResid2 .  ^2) ; 

SStotal2=(length(Y2)-l)*var(Y2); 

rsq2=l-SSresid2./SStotal2; 

adj_rsq2=  1  -S  Sresid2  ./S  Stotal2  *(length(  Y  2)- 1 )  ./(length(Y2)- . . . 
length(ReeF  req2)) ; 

%See.  5.3.  Baseline  Curve  2.5  deg  Below  Beam  Center 
lat3=(lat-2.5)*pi/180;  %deg  to  rad 
Nphi3=Re/sqrt(  1  -e.^2  *  (sin(lat3))  .^2); 
JamPosX3=a-Nphi3*cos(lat3)*cos(  Angle); 
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JamPosY3=Nphi2*cos(lat3)*sin(  Angle); 

JamPosZ3=(  1  -e  A2)*Nphi3  *sin(lat3); 

JamPos3=sqrt(  JamPosX3  A2+JamPosY  3  A2+JamPosZ3  A2); 

Rdot3=(diff(JamPos3))/dt; 

Y  3 =(Rdot3/c)  *  Fnot  J  am+F  not  Jam; 
plot  (Tbase,  Y3); 

hold  on 

RecFreq3=polyfit(Tbase,Y3,3); 

XTrack3=polyval(RecFreq3,maxTime); 

YJamFit3=polyval(RecFreq3, Tbase); 

Y  JaniResid3=Y3  -  Y  JamFit3 ; 

S  Sresid3 =sum(  Y  J  amResid3  A2) ; 

SStotal3=(length(Y3)-l)*var(Y3); 

rsq3=l-SSresid3./SStotal3; 

adj_rsq3= 1  -S  Sresid3  ./S  StotaB  *(length(  Y3)- 1 )  ./(length(Y  3 )- . . . 
length(ReeF  req3 )) ; 

5lc  sic  sic  sic  *  ![C  !lc  sic  *  ![C  !lc  ![C  *  #  !lc  ![C  *  sic  sic  sic  *  ![C  !lc  sic  *  sic  sic  ![C  *  sic  !lc  sic  >1C  ^  sic  sic  *  ![C  !lc  sic  *  sic  !lc  sic  !lc  ![C  !lc  sic  !lc  sic  !lc  sic  *  ![C  !lc  sic  *  sic  !lc  sic  *  ![C  !lc  sic  !lc  sic  !lc  sic  *  ![C 

%See.  5.4.  Baseline  Curve  2  deg  Below  Beam  Center 
lat4=(lat-2)*pi/180;  %deg  to  rad 
Nphi4=Re/sqrt(  1  -e.^2  *  (sm(lat4))  .^2); 

JamPosX4=a-Nphi4*cos(lat4)*cos(  Angle); 

JamPosY4=Nplii2*eos(lat4)*sm(  Angle); 

JamPosZ4=(l-e.^2)*Nphi4*sin(lat4); 

JamPos4=sqrt(JamPosX4.^2+JamPosY4.^2+JamPosZ4.^2); 

Rdot4=(diff(JamPos4))/dt; 

Y  4=(Rdot4/ e)  *  Fnot  J  am+F  not  Jam; 
plot  (Tbase,  Y4); 

hold  on 

ReoFreq4=polyfit(Tbase,Y4,3); 

XTrack4=polyval(ReoFreq4,maxTime); 

YJamFit4=polyval(ReoFreq4, Tbase); 

YJamResid4=Y4-YJamFit4; 

S  Sresid4=sum(  Y  J  amResid4 .  ^2) ; 

S  S  total4=(length(  Y  4)-l)*var(Y4); 
rsq4=l-SSresid4./SStotal4; 

adj_rsq4=  1  -S  Sresid4  ./S  StotaM  *(length(  Y  4)- 1 )  ./(length(Y  4)- . . . 
length(ReeF  req4)) ; 

%See.  5.5.  Baseline  Curve  1.5  deg  Below  Beam  Center 
lat5=(lat-1.5)*pi/180;  %deg  to  rad 
Nphi5=Re/sqrt(  1  -e.^2  *  (sm(lat5))  .^2); 

JamPosX5=a-Nphi5*oos(lat5)*oos(  Angle); 

JamPosY5=Nphi2*cos(lat5)*sm(  Angle); 
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JamPosZ5=(  1  -e.^2)*Nphi5  *sin(lat5); 

JamPos5=sqrt(  JamPosX5  Zl+JamPosY  5  Z2+JamPosZ5  .^2); 
Rdot5=(diff(JamPos5))/dt; 

Y  5 =(Rdot5/c)  *  Fnot  J  am+F  not  Jam; 
plot  (Tbase,  Y5); 

hold  on 

RecFreq5=polyfit(Tbase,Y5,3); 
XTrack5=polyval(RecFreq5,maxTime); 
YJamFit5=polyval(RecFreq5, Tbase); 

Y  JamResid5=Y5  -  Y  JamFitS ; 

S  SresidS =sum(  Y  J  amResidS .  ^2) ; 

SStotal5=(length(Y  5)- 1  )*  var(Y  5); 
rsq5=l-SSresid5./SStotal5; 

adj_rsq5= 1  -S  SresidS  ./S  Stotal5  *(length(  Y  5)- 1 )  ./(length(Y  5  )- . . . 
length(RecF  req5  )) ; 

%Sec.  5.6.  Baseline  Curve  1  deg  Below  Beam  Center 
lat6=(lat-l)*pi/180;  %deg  to  rad 
Nphi6=Re/sqrt(  1  -e.^2  *  (sm(lat6))  .^2); 
JamPosX6=a-Nphi6*cos(lat6)*cos(  Angle); 
JamPosY6=Nphi2*eos(lat6)*sm(  Angle); 
JamPosZ6=(l-e.^2)*Nphi6*sin(lat6); 
JamPos6=sqrt(JamPosX6.^2+JamPosY6.^2+JamPosZ6.^2); 
Rdot6=(diff(JamPos6))/dt; 

Y  6=(Rdot6/ c)  *  Fnot  J  am+F  not  Jam; 
plot  (Tbase,  Y6); 

hold  on 

ReoFreq6=polyfit(Tbase,Y6,3); 
XTraok6=polyval(ReoFreq6,maxTime); 
YJamFit6=polyval(ReoFreq6, Tbase); 
YJamResid6=Y6-YJamFit6; 

S  Sresid6=sum(  Y  J  amResidb .  ^2) ; 

SStotal6=(length(Y  6)- 1  )*  var(Y  6); 
rsq6=l-SSresid6./SStotal6; 

adj_rsq6=  1  -S  Sresidb  ./S  Stotalb  *(length(  Y  6)- 1 )  ./(length(Y  6)- . . . 
length(RecF  req6)) ; 

%See.  5.7.  Baseline  Curve  .5  deg  Below  Beam  Center 
lat7=(lat-.5)*pi/180;  %deg  to  rad 
Nphi7=Re/sqrt(  1  -e.^2  *  (sm(lat7))  .^2); 
JamPosX7=a-Nphi7*cos(lat7)*cos(  Angle); 
JamPosY7=Nphi2*cos(lat7)*sm(  Angle); 
JamPosZ7=(l-e.^2)*Nphi7*sin(lat7); 
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JamPos7=sqrt(JamPosX7.^2+JamPosY7.^2+JamPosZ7.^2); 

Rdot7=(diff(JamPos7))/dt; 

Y  7=(Rdot7/c)  *  Fnot  J  am+F  not  J  am; 
plot  (Tbase,  Y7); 

hold  on 

RecFreq7=polyfit(Tbase,Y7,3); 
XTrack7=polyval(RecFreq7,maxTime); 
YJaniFit7=polyval(RecFreq7, Tbase); 
YJamResid7=Y7-YJamFit7; 

S  Sresid7=sum(  Y  J  amResid7 .  ^2) ; 

SStotal7=(length(Y7)-l)*var(Y7); 

rsq7=l-SSresid7./SStotal7; 

adj_rsq7=  1  -S  Sresid7  ./S  Stotal7  *(length(  Y7)- 1 )  ./(length(Y7)- . . . 
length(RecF  req7)) ; 

%Sec.  5.8.  Baseline  Curve  Beam  Center 
lat8=(lat)*pi/180;  %deg  to  rad 
Nphi8=Re^qrt(  1  -e.^2  *  (sm(lat8))  .^2); 
JamPosX8=a-Nphi8*oos(lat8)*oos(  Angle); 
JamPosY8=Nplii2*eos(lat8)*sm(  Angle); 
JamPosZ8=(l-e.^2)*Nphi8*sin(lat8); 
JamPos8=sqrt(JamPosX8.^2+JamPosY8.^2+JamPosZ8.^2); 
Rdot8=(diff(JamPos8))/dt; 

Y  8=(Rdot8/ e)  *  Fnot  J  am+F  not  J  am; 
plot  (Tbase,  Y8); 

hold  on 

ReoFreq8=polyfit(Tbase,Y8,3); 
XTraok8=polyval(ReoFreq8,maxTime); 
YJamFit8=polyval(ReoFreq8, Tbase); 
YJamResid8=Y8-YJamFit8; 

SSresid8=sum(YJamResid8.^2); 

SStotal8=(length(Y8)-l)*var(Y8); 

rsq8=l-SSresid8./SStotal8; 

adj_rsq8= 1  -S  Sresid8  ./S  Stotal8  *(length(  Y  8)- 1 )  ./(length(Y8)- . . . 
length(ReeFreq8)); 

%See.  5.9.  Baseline  Curve  .5  deg  Above  Beam  Center 
lat9=(lat+.5)*pi/180;  %deg  to  rad 
Nphi9=Re/sqrt(  1  -e.^2  *  (sm(lat9))  .^2); 
JamPosX9=a-Nphi9*oos(lat9)*oos(  Angle); 
JamPosY9=Nphi2*eos(lat9)*sm(  Angle); 
JamPosZ9=(l-e.^2)*Nphi9*sin(lat9); 
JamPos9=sqrt(JamPosX9.^2+JamPosY9.^2+JamPosZ9.^2); 
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Rdot9=(diff(JamPos9))/dt; 

Y  9=(Rdot9/ c)  *  Fnot  J  am+F  not  J  am; 
plot  (Tbase,  Y9); 

hold  on 

RecFreq9=polyfit(Tbase,Y9,3); 
XTrack9=polyval(RecFreq9,maxTime); 
YJaniFit9=polyval(RecFreq9, Tbase); 

Y  JaniResid9=Y  9- Y  JamFit9; 

S  Sresid9=sum(  Y  J  amResid9 .  ^2) ; 

S  S  total9=(length(  Y  9)-l)*var(Y9); 
rsq9=l-SSresid9./SStotal9; 

adj_rsq9=  1  -S  Sresid9  ./S  Stotal9  *(length(  Y  9)- 1 )  ./(length(Y  9)-... 
length(RecF  req9)) ; 

%Sec.  5.10.  Baseline  Curve  1  deg  Above  Beam  Center 
latl0=(lat+l)*pi/180;  %deg  to  rad 
Nphi  1 0=Re/sqrt(  l-e.^2*(sin(latl0)).^2); 

JamPosX  1 0=a-Nphi  1 0  *  cos(lat  1 0)  *  cos(  Angle); 
JamPosY10=Nphi2*eos(latl0)*sm(  Angle); 
JamPosZ10=(l-e.^2)*Nphil0*sm(latl0); 
JamPosl0=sqrt(JamPosX10.^2+JamPosY10.^2+JamPosZ10.^2); 
RdotlO=(diff(JamPoslO))/dt; 

Y 1 0=(Rdot  1 0/e)  *  Fnot  J  am+F  not  Jam; 
plot  (Tbase,  YIO); 
hold  on 

ReeFreql  0=polyfit(Tbase,  Y  10,3); 
XTracklO=polyval(RecFreqlO,maxTime); 

Y  JamFit  1 0=poly  val(RecFreq  1 0, Tbase); 

YJamResidl  0=Y  1 0- YJamFit  1 0; 

SSresid  1 0=sum(  Y  J  amResid  1 0  .^2); 

SStotal  1 0=(length(  Y 1 0)- 1 )  *  var(Y  1 0); 
rsql0=l-SSresidl0./SStotall0; 

adj_rsql0=l-SSresidl0./SStotall0*(length(Y10)-l)./(length(Y10)-... 
length(ReeF  req  1 0)) ; 

%See.  5.11.  Baseline  Curve  1.5  deg  Above  Beam  Center 
latll=(lat+1.5)*pi/180;  %deg  to  rad 
Nphil  l=Re/sqrt(l-e.^2*(sm(latl  1)).^2); 

JamPosXl  l=a-Nphil  1  *oos(latl  l)*cos(  Angle); 

JamPosY  1  l=Nphi2*cos(latl  l)*sm(Angle); 

JamPosZl  l=(l-e.^2)*Nphil  l*sm(latl  1); 

JamPosl  l=sqrt(JamPosXl  1  .^2+JamPosYl  1  .^2+JamPosZl  1  .^2); 
Rdotl  l=(diff(JamPosl  l))/dt; 


77 


Y 1  l=(Rdotl  l/c)*FnotJam+FnotJam; 
plot  (Tbase,  Yll); 
hold  on 

RecFreql  l=polyfit(Tbase,Yl  1,3); 

XTrackl  l=polyval(RecFreql  l,maxTime); 

YJamFitl  l=polyval(RecFreql  1, Tbase); 

YJamResidl  1=Y1 1 -YJamFitl  1; 

SSresidl  l=sum(  YJamResidl  1  .^2); 

SStotall  l=(length(Yl  l)-l)*var(Yl  1); 
rsql  l=l-SSresidl  1  ./SStotall  1; 

adj_rsql  l=l-SSresidl  1. /SStotall  l*(length(Yl  l)-l)./(length(Yl  1)-... 
length(ReoFreql  1)); 

%Seo.  5.12.  Baseline  Curve  2  deg  Above  Beam  Center 
latl2=(lat+2)*pi/180;  %deg  to  rad 
Nphi  1 2=Re/sqrt(  l-e.^2*(sin(latl2)).^2); 

JamPosX  1 2=a-Nphi  1 2  *  cos(lat  1 2)  *  cos(  Angle); 
JamPosY12=Nphi2*eos(latl2)*sm(Angle); 

JamPosZ  1 2=(  1  -e.^2)*Nphi  1 2*sm(lat  1 2); 

JamPosl2=sqrt(JamPosX12.^2+JamPosY12.^2+JamPosZ12.^2); 

Rdotl2=(diff(JamPosl2))/dt; 

Y12=(Rdotl2/e)*FnotJam+FnotJam; 
plot  (Tbase,  Y12); 
hold  on 

ReoFreql2=polyfit(Tbase,Y12,3); 

XTraokl2=polyval(ReoFreql2,maxTime); 

Y  JamFit  1 2=poly  val(ReoFreq  1 2, Tbase); 

Y  JamResid  1 2=Y  1 2- Y  JamFit  12; 

SSresid  1 2=sum(  Y  J  amResid  1 2  .^2); 

SStotal  1 2=(length(Y  1 2)- 1  )*  var(Y  1 2); 
rsql2=l-SSresidl2./SStotall2; 

adj_rsql2=l-SSresidl2./SStotall2*(length(Y12)-l)./(length(Y12)-... 
length(RecF  req  1 2)) ; 

%See.  5.13.  Baseline  Curve  2.5  deg  Above  Beam  Center 
latl3=(lat+2.5)*pi/180;  %deg  to  rad 
Nphil3=Re/sqrt(l-e.^2*(sm(latl3)).^2); 
JamPosX13=a-Nphil3*cos(latl3)*cos(Angle); 
JamPosY13=Nphi2*cos(latl3)*sm(  Angle); 

JamPosZ  1 3=(  1  -e.^2)*Nphi  1 3  *sm(lat  1 3); 

JamPos  1 3=sqrt(  JamPosXl  3  .^2+JamPosY  1 3  .^2+JamPosZ  1 3  .^2); 
Rdotl3=(diff(JamPosl3))/dt; 

Y13=(Rdotl3/c)*FnotJam+FnotJam; 
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plot  (Tbase,  Y13); 
hold  on 

RecFreql  3=polyfit(Tbase,  Y  13,3); 

XTrackl3=polyval(RecFreql3,maxTime); 

Y  JamFit  1 3=poly  val(RecFreq  1 3 , Tbase); 

Y  JamResid  1 3=Y  1 3  -  Y  JamFit  1 3 ; 

SSresid  1 3=sum(  Y  J  amResid  1 3  .^2); 

SStotall3=(length(Y13)-l)*var(Y13); 

rsql3=l-SSresidl3./SStotall3; 

adj_rsql3=l-SSresidl3./SStotall3*(length(Y13)-l)./(length(Y13)-... 
length(RecF  req  1 3 )) ; 

=1=  *  =1=  He  *  =[c  =lc  He  =1=  =[c  =1=  He  *  =[c  =lc  He  =1=  He  =1=  He  *  =[c  =lc  He  =1=  =[c  =1=  =[c  >1=  #  =lc  =[c  =1=  He  =1=  He  *  =[c  =lc  He  =1=  *  =1=  He  =1=  =[c  =lc  He  =1=  He  =1=  He  *  =[c  =lc  He  *  He  =1=  He  *  =[c  =lc  He  =lc  He  =1=  He  *  =[c 

%Sec.  5.14.  Baseline  Curve  3  deg  Above  Beam  Center 
latl4=(lat+3)*pi/180;  %deg  to  rad 
Nphi  1 4=Re/sqrt(  l-e.^2*(sin(latl4)).^2); 

JamPosX14=a-Nphil4*oos(latl4)*cos(  Angle); 
JamPosY14=Nphi2*cos(latl4)*sin(Angle); 

JamPosZ  1 4=(  1  -e.^2)*Nphi  1 4*sin(lat  1 4); 

JamPosl4=sqrt(JamPosX14.^2+JamPosY14.^2+JamPosZ14.^2); 

Rdotl4=(diff(JamPosl4))/dt; 

Y14=(Rdotl4/c)*FnotJam+FnotJam; 
plot  (Tbase,  Y14); 
hold  on 

ReoFreql4=polyfit(Tbase,Y14,3); 

XTrackl4=polyval(RecFreql4,maxTime); 

Y  JamFit  1 4=poly  val(ReoFreq  1 4, Tbase); 

YJamResidl4=Y14-YJamFitl4; 

SSresid  1 4=sum(  Y  J  amResid  1 4  .^2); 

SStotal  1 4=(length(  Y 1 4)- 1 )  *  var(Y  1 4); 
rsql4=l-SSresidl4./SStotall4; 

adj_rsq  1 4=  1  -SSresid  1 4  ./S  Stotal  1 4  *(length(  Y 1 4)- 1 )  ./(length(Y  14)-... 
length(ReeF  req  1 4)) ; 

%Sec.  5.15.  Baseline  Curve  3.5  deg  Above  Beam  Center 
latl5=(lat+3.5)*pi/180;  %deg  to  rad 
Nphil5=Re/sqrt(l-e.^2*(sm(latl5)).^2); 

JamPosX15=a-Nphil5*cos(latl5)*cos(Angle); 

JamPosY15=Nphi2*oos(latl5)*sm(  Angle); 

JamPosZ  1 5=(  1  -e.^2)*Nphi  1 5  *sm(lat  1 5); 

JamPos  1 5=sqrt(  JamPosXl  5  .^2+JamPosY  1 5  .^2+JamPosZ  1 5  .^2); 
Rdotl5=(diff(JamPosl5))/dt; 

Y15=(Rdotl5/o)*FnotJam+FnotJam; 
plot  (Tbase,  Y15); 
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hold  on 

RecFreql  5=polyfit(Tbase,  Y  15,3); 

XTrackl5=polyval(RecFreql5,maxTime); 

Y  JamFit  1 5=poly  val(RecFreq  1 5  ,Tbase); 

YJamResidl  5=Y1 5-YJamFitl  5; 

SSresid  1 5=sum(  Y  J  amResid  1 5  .^2); 

SStotall5=(length(Y15)-l)*var(Y15); 

rsql5=l-SSresidl5./SStotall5; 

adj_rsql5=l-SSresidl5./SStotall5*(length(Y15)-l)./(length(Y15)-... 

length(RecFreql5)); 

legend  Jam  -3.5deg  -3deg  -2.5deg  -2deg  -1.5deg  -Ideg  -.5deg  BeamCenter... 

+.5deg  +ldeg  +1.5deg  +2deg  +2.5deg  +3deg  +3.5deg; 
hold  off 

%Seo.  6.  Interpolation  of  Jammer  Latitude  from  Base  Curves 

5lc  sic  sic  sic  5lc  ![C  sic  sic  *  ^  >1C  sic  sic  #  !lc  sic  !lc  ![C  !lc  ![C  5lc  #  sic  sic  *  ![C  !lc  ![C  *  ![C  !lc  sic  !lc  sic  !lc  sic  !lc  ^  !lc  sic  *  ![C  sic  ![C  5lc  ![C  ![C  !lc  sic  sic  sic  sic  #  !lc  sic  !lc  sic  !lc  ![C  *  ![C  !lc  sic  !lc  sic  !lc  sic  *  ![C 

YT=[latl  lat2  lat3  lat4  lat5  lat6  latV  lat8  lat9  latlO  latll  latl2  latl3... 

Iatl4  latl5]; 

XT=[XTraekl  XTraek2  XTraek3  XTraek4  XTraek5  XTrackb  XTraekV  XTrack8... 

XTraek9  XTraeklO  XTraekl  1  XTraekl2  XTraekl3  XTraekl4  XTraekl5]; 
XTSolve=polyfit(XT,YT,4); 

XT  rkF  it=poly  val(XT  S  olve,XT) ; 

XTrkResid=YT-XTrkFit; 

S  SresidXT=sum(XT  rkResid.  ^2) ; 

SStotalXT=(length(YT)- 1  )*  var(YT); 
rsqXT=l-SSresidXT./SStotalXT; 

adj_rsqXT=  1  -S  SresidXT  ./S  StotalXT*  (length(YT)- 1 )  ./(length(  YT)- . . . 

length(XTSolve)); 

%figure 

%plot  (XT,  YT); 
hold  off 

LatRad=polyval(XTSolve,XTraokJam);  %Latitude  in  radians 
LatEst=LatRad*  180/pi  %radians  to  degrees 

%See.  7.  Display  Estimated  Position  on  Map  (Red  X) 
figure 

axesm  eekert4; 
framem;  gridm; 
axis  on 

geoshow('landareas.shp','Eaoeoolor','y'); 

hnem(EatEst,EongEst,'EineStyle','none','LineWidth',2,'Color','r',.-. 
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'Marker', 'x','MarkerSize',  1 0); 
hold  on 

%Sec.  8.  Calculate  Error  from  Actual  Position 

[LatDist]=distance(ActLat,  ActLong,  LatEst,  ActEong,  ellipsoid)/! 000  %km 
[LongDist]=distance(ActEat,  ActEong, ActEat,  EongEst,  ellipsoid)/! 000  %km 
[Dist]=distance(ActEat,  ActEong,  EatEst,  EongEst,  e!!ipsoid)/!000  %km 

%Sec.  9.  Eongitude  Error  Correction/P!ot  Based  on  Obtained  resuits  (Red  +) 

=1=  He  =1=  He  *  =[c  He  =1=  He  =1=  He  =1=  #  =lc  =[c  =1=  He  =1=  =[c  *  He  He  He  =1=  =[c  =lc  He  >1=  =lc  =lc  He  =1=  *  =1=  He  >1=  =[c  He  =1=  *  =1=  He  =1=  =[c  =lc  He  =1=  He  =1=  He  =1=  =lc  He  =1=  =[c  =1=  =[c  *  =[c  =lc  =[c  =1=  He  =1=  He  >1=  =[c 

EWErr=[-!.485  -2.838  -5.900  -10.931]; 

!ats=[15  30  45  60]; 

EWCorfit=po!yfit(!ats,  EWErr,  2); 

EWCor=po!yva!(EWCorfit,  EatEst); 

EongCor=Sat!nitEong-l-(RotV*(Sat!nitTime+AbeamTime)*180/pi)-(EWCor*RotV*... 

180/pi)+.0649 

[LongCorDist]=distance(ActEat,  ActEong,  ActEat,  EongCor,  ellipsoid)/! 000  %km 
[CorDist]=distance( ActEat,  ActEong,  EatEst,  EongCor,  ellipsoid)/! 000  %km 
bnem(EatEst,EongCor,'EineStyle','none','EineWidth',2,'Color','r',... 

'Marker','-!-', 'MarkerSize',  1 0); 
hold  on 

%Sec.  10.  Plot  Actual  Eocation  on  Map  for  Comparison  (Green  X) 

bnem(  ActEat,  ActEong,'EineStyle','none','EineWidth',2,'Color','g',... 

'Marker','X', 'MarkerSize',  1 0); 
hold  off 
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